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Abstract: Various holographic approaches to QCD in five dimensions are explored 
using input both from the putative five-dimensional non-critical string theory as well 
as QCD. It is argued that a gravity theory in five dimensions coupled to a dilaton and 
an axion may capture the important qualitative features of pure YM theory. A part of 
the effects of higher a'-corrections is resummed into a dilaton potential. The potential 
is shown to be in one-to-one correspondence with the exact /5-function of QCD, and 
its knowledge determines the full structure of the vacuum solution. The geometry 
near the UV boundary is that of AdSs with logarithmic corrections reflecting the 
asymptotic freedom of QCD. We find that all relevant confining backgrounds have an 
IR singularity of the "good" kind that allows unambiguous spectrum computations. 
Near the singularity the 't Hooft coupling is driven to infinity. Asymptotically linear 
glueball masses can also be achieved. The classification of all confining asymptotics, 
the associated glueball spectra and meson dynamics are addressed in a companion 
paper |ArXiv:0707.1349| 
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1. Introduction and conclusions 

The correspondence between large- N c gauge theory and string theory first advocated 
by 't Hooft [Q has taken a novel (and mostly unexpected) form after the Maldacena 
conjecture, 0. Although the most precise version of the correspondence was formu- 
lated for the (conformally invariant) maximally supersymmetric Yang Mills theory 
in four dimensions, in several other cases, holographic techniques proved useful, and 
several gravitational duals describing confining theories in the IR have been proposed. 
One concrete application of holographic duality is to establish a novel quantitative 
understanding of QCD, in particular of phenomena where strong IR physics is impor- 
tant. These include confinement, chiral symmetry breaking, as well as quantitative 
issues about the meson and baryon spectrum and their interactions. 

Critical holographic models obtained as solutions of ten-dimensional string the- 
ory H ||, H have been successful in their qualitative description of confinement and 
other IR dynamics, including a semi- quantitative agreement of the glueball spectrum 
[0] and thermal properties [[|. At the same time, the theories contain Kaluza-Klein 
modes, not expected in QCD, with KK masses of the same order as the dynamical 
scale of the gauge theory. Above this scale the theories deviate from QCD. In the 
solution associated to non-supersymmetric wrapped D4 branes, || the UV comple- 
tion should be thought of as a 6D theory on M5 branes. In the Chamsedinne-Volkov 
solution[|J, interpreted holographically by Maldacena and Nunez, [[|, the theory be- 
comes the 6d theory obtained by wrapping NS5 branes on a two-sphere. Finally in 
the Klebanov-Strassler solution, || the theory passes through a large set of Seiberg 
dualities to end up with a non-conformal quiver theory in the UV. 
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An obvious way to avoid the extra KK modes is to consider non-critical string 
theories. Indeed, if our general intuition about holography is correct, QCD should 
be described by a string theory living in 5 dimensions. However, it is expected from 
general arguments that non-critical string duals of large- N c gauge theories will have 
background curvatures and other invariants of the same order as the string scale 
and an a 1 expansion is not expected to be a priori a reliable approximation. Such 
an expectation is partly based on the fact that the part of the potential due to the 
central charge deficit is of order 0(1) in t~ 2 units. One may also give a more rigorous 
statement || for supersymmetric theories, based on the conformal anomaly @. 

Despite the hostile environment of non-critical theory, several attempts have 
been made to understand holographic physics in this regime, based on two-derivative 
gravitational actions, [JTOj, [11], §, ^2| . Moreover, such attempts were seconded by 
studies of the physics of probe branes in non-critical backgrounds [13|, [14], and non- 
critical orientifolds that provide side information on the holographic approaches. 
Most of these attempts focused on 4d CFTs with M = 1 supersymmetry [[□], Q 
or without it [K], |8|, [HJ. The rational has been, that although a' corrections are 
generically expected to be substantial, qualitative information should be (mostly) 
reliable. Moreover, the high symmetry may guarantee that some quantities could 
be reliably computed. Part of this intuition stems from exact solutions to 2d WZW 
CFTs where ratios of conformal weights can be calculated classically, and are exact 
to all orders in the a' expansion. 

A different and more phenomenological approach was in the meantime developed, 
and is now known as AdS/QCD. The original idea was formulated in W7\ and it 
was successfully applied to the meson sector in |TB], |HJ. The bulk gravitational 
background consists of a slice of AdS^, and a constant dilaton. There is a UV and 
an IR cutoff. Moreover, the confining IR physics is imposed by boundary conditions 
at the IR boundary. This approach seems very crude and indeed it is, when applied 
to the pure glue sector. However it has been partly successful in studying meson 
physics despite the fact that the dynamics driving chiral symmetry breaking must 
be imposed by hand via IR boundary conditions. 

In this paper and its companion p0| , we will investigate a mostly phenomenolog- 
ical approach that runs somewhere in-between non-critical string theory approaches 
and AdS/QCD. On the one hand we will investigate and motivate what kind of effec- 
tive theory we expect to describe QCD, based on our understanding of string theory. 
On the other hand we would like to match this with what we expect from QCD in 
the UV and the IR. 

There is intriguing evidence from studies of QCD, that high-dimension operators, 
that should be associated to stringy modes, are not very important for (some of the) 
physics at short distances. This is suggested by the success of SVZ sum rules, [^TJ 
which tie together the UV physics with the IR physics in QCD. A counter-argument 
relies on the fact that as we understand from critical holography, the large 't Hooft 
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coupling limit A ^> 1 is necessary to suppress higher a' corrections. Therefore, in 
QCD, as this coupling is driven to zero in the UV, such corrections will become 
dominant. However, as we argue below, progress can be achieved in this direction 
despite the difficulties. 

Pure 4d YM at large N c is expected to be dual to a string theory living in 
5 dimensions. The relevant low-lying fields are expected to be dual to the lowest 
dimension operators, namely the graviton (dual to T^ u ~ Tr[F 2 v — ^S^F 2 }), the 
dilaton (dual to Tr[F 2 }) and the RR axion a (dual to Tr[F A F]). Moreover, in 
the string frame, the theory has a very simple dilaton potential ~ 5c that reflects 
the fact that the associated string theory is non-critical. The theory however must 
contain a RR four-form whose flux seeds the D% branes and therefore generates the 
U(N C ) gauge group. The presence of such a form is compatible with the spectrum 
of operators we advocate, as in 5 dimensions a four-form is non-propagating. By 
integrating out the four-form we generate new terms in the dilaton potential. This 
is shown in section |^ and appendix [TJ. Therefore, taking into account higher a' 
terms involving the four-form only, generates a non-trivial potential for the dilaton 
(j). Hence, the special nature of the four-form implies that its field strength, although 
it carries powers of a', effectively carries no derivatives, and it should therefore be 
incorporated in the potential. 

There are two properties that this potential satisfies. The first is that it is a 
nontrivial function of the 't Hooft coupling A ~ N c e^, (A thus defined, is expected to 
remain finite in the large- N c limit). The second is that (a) it has a regular expansion 
around A = (b) it vanishes as At at A = (in the Einstein frame). To next to 
leading order we expect 

V(\) = \l (5c-\ 2 + ---) (1.1) 

This potential was established in || and is generic in dimensions 5, 4 and 3. The 
fact that it vanishes at A = implies that the asymptotic geometry in the UV does 
not approach AdSs, as we might expect. This is analyzed in detail in appendix O. 
However the potential above has another intriguing property: it has an AdS minimum 
at a finite value of the 't Hooft coupling, and therefore there is a related AdSs solution 
with fixed finite 't Hooft coupling. The expectation that this AdSs may be the UV 
limit of a non-trivial solution with (logarithmically) running A seems however to fail 
as explained in appendix 0. Therefore the only "correct" weak-coupling asymptotics 
of the dilaton potential must satisfy 

lim V (A) = V ^ (1.2) 

This guarantees the existence of a solution which is asymptotically AdSs near the 
UV boundary, provided the 't Hooft coupling vanishes there. Moreover, if we wish 
that in the same UV region, the 't Hooft coupling runs logarithmically with energy, 
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then the weak coupling expansion of the potential must be of the form 

V(X) = V f l + f>„(AT) (1-3) 



n=l 



with a some positive real number that can be shifted to the wave function renor- 
malization of the dilaton. As suggested by the perturbative QCD /3-function we 
will select a = 1. As we show in section |2| and appendix [A], this expansion of the 
potential is equivalent to the perturbative expansion of the QCD /3-function and the 
coefficients V n can be directly related to the perturbative /3-function coefficients. 
There are two obvious questions that accompany the discussion above: 

1. What is the origin of the non-zero constant Vq? 

2. In the arguments above we have neglected higher a'-corrections involving the 
curvature and derivatives of the scalars. 



As we indicate in appendix |B.3| , a non-zero constant Vq in the effective dilaton 
potential may be generated from the higher-curvature corrections. Although this also 
implies that higher cV-corrections cannot be neglected after all, we would like to take 
the bold step and assume that the only quantitatively important part of this class of 
a 1 corrections is just to generate a non-zero Vq as well as the rest of the terms of the 
dilaton potential. This step, has as a consequence that we will eventually deal with 
a two-derivative effective action but with a general dilaton potential satisfying ( |1.3| ) 
at weak 't Hooft coupling. Although this approach cannot be rigorously defended, 
our attitude is exploratory. An extra motivation comes from the success of the SVZ 
sum rules [ 21[1 . Moreover, as we show below, our results are encouraging. 



To summarize the discussion above, our starting point is an action (in the Ein- 
stein frame) of the form 



1.4) 



S = M 3 N* J d 5 x^j 

with the potential V(X) having around A = an asymptotic expansion of the form 



H). 

The axion kinetic term in (|1.4j ) deserves some comments. In a way similar with 
the potential, we expect a non-trivial function of the dilaton Z(X) multiplying its 
kinetic term. 1 Moreover, as explained in detail in appendix |B.1| , the axion kinetic 
term is of order 0(1/N^) compared with the rest of the terms in ( |1.4j ). This is 
due in string theory to the fact that the axion is a RR field and therefore has a 
suppressed dilaton dependence. 2 The same argument indicates that terms involving 



1 A priori, similar functions might also multiply the Einstein term and the dilaton kinetic term. 
Such terms can however be removed by Weyl-rescaling of the metric and redefining A. 

2 This is not valid for the RR four-form as its field strength is linear in the number of colors N c . 
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higher powers of (da) 2 will be further suppressed at large N c . Therefore, the a'- 
expansion of the axion terms is effectively an 1/N C expansion. The perturbative 
asymptotics of Z(X) therefore should be 

Z(X) = Z a + 0(X 2 ) , A-^0 (1.5) 

String theory gives a A 2 contribution to leading order. However as in the case of 
the potential, higher derivative corrections are expected to generate also a constant 
piece. This turns out to be in agreement with perturbative QCD. 

As shown in section 4 of [201, Z(X) determines an analogue of a /3-function for 



the QCD ^-parameter. This interpretation has however caveats that are discussed in 
the same section. It is obvious that while searching for the solution of the equations 
of motion stemming from the action (|1.4j) , that will describe the QCD vacuum, the 
axion contribution can be neglected to leading order in 1/N C . Once this solution is 
found, the axion equation of motion can be solved in order to determine the profile 
of the QCD ^-parameter and its associated physics. 

Therefore, to leading order in 1/N C , the vacuum structure of QCD, is captured 
by a solution of ( |1.4j ) with the appropriate AdS<$ asymptotics. All the properties 
of the 4D gauge theory depend on a single function of the 't Hooft coupling, the 
superpotential W(X), defined as 3 : 

i 2 fdW\ 2 , 64 2 



This dependence can be mapped (up to the overall AdS length £, related to V by 
= |i) to the (exact) /3-function of QCD, /3(A). In particular, for small A, all 
coefficients of the superpotential can be mapped to the perturbative terms of the 
/3-function. Concretely, the following relation holds: 

= _9 a rflog^A) 
4 d\ 



Moreover fitting to QCD data we learn that, 

£-6.26 4 (1.8) 

We expect three integration constants for the equations of motion of the metric 
and A. One of them fixes W(X) as a solution of eq. ( |1.6|) . We show in [20 that this 



integration constant is completely fixed by using asymptotic freedom as an input 
from the gauge theory. Among the second and third integration constants that arise 



3 Generically, passing from the potential to the superpotential and the associated first order 
equations involves an arbitrary constant that can be thought of as the single constant of integration 



of (1.6). However in our case it turns out that the relevant solution is special, and does not have 



this extra dependence on an arbitrary constant. This is explained in appendix E of |2C 
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from the first-order differential equations for A and A } only a single combination 
remains due to the reparametrization invariance. This remaining integration constant 



amounts to a definition of Aq C d- ^From eq. ( |L7| ) it is apparent that the exact (3- 
function determines completely the geometry up to a choice of Aq C d- 

One result, in our two-derivative approach, concerns the investigation of potential 
confining backgrounds in the IR. We choose the conformal coordinate system and 
write the Einstein metric as 

ds 2 = e 2A{r \dr 2 + rj^dx^dx 11 ) (1.9) 

with the AdSs boundary at r = 0. As we show, e j4(r - ) monotonically decreases from 
oo at the UV boundary, to 0, or to a finite value, in the IR (see section 3.2). We 
therefore take the scale factor in the Einstein frame as our definition for the energy 
of the gauge theory. 

E = e A{r) (1.10) 

In the case when the Einstein frame scale factor remains finite at the IR singularity, 
this indicates that the dual theory is defined only above a certain energy. Such 
models however fail to satisfy some of the properties that are believed to hold in 
confining theories, namely the screening of the magnetic color charges. 

In our analysis of the non-perturbative regime, we give a general classification of 
possible IR geometries according to their confining properties. We charac- 
terization of confinement, the Wilson loop area law: by "confining," we label those 
backgrounds for which the QCD string (identified with the fundamental string) has 
a finite tension. We then analyze various properties of confining backgrounds. Here 
is a summary of our findings: 

• We provide a relation between the /5-function of the gauge theory defined by 
an infinite series expansion in the 't Hooft coupling, and the scalar potential 
(or rather the superpotential) of the gravitation theory defined by a similar 
expansion in the dilaton. 

• We study certain a' corrections to the scalar potential. In particular, we find 
that in the identification of the /^-function and scalar potential coefficients, 
the first two scheme-independent coefficients of the /3-function receive no a' 
corrections. Moreover we find that the a' corrections to the identification of the 
energy scale with the scale factor of the metric can be set to zero in a particular 
scheme for computing the higher order /3-function coefficients, 62, 63, ... . 

• We show that all confining backgrounds have a singularity in the Einstein 
metric at some value r$ of the r coordinate. There are two distinct cases for 
the position r of the IR singularity. One possibility is finite r . The other 
is r = 00. The IR singularity is always of the "good kind" ||22|| : fluctuation 
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spectra of various fields are well defined and are not affected directly by the 
presence of the singularity. 



For regular dilaton potentials 4 , the 't Hooft coupling A always becomes infinite 
at the IR singularity. 

In the case r = oo, the string frame metric is not only regular at the IR 
singularity but its curvature also vanishes. Put otherwise, in the string frame, 
the IR singularity is only due to the diverging 't Hooft coupling constant. This 
suggests that the supergravity approximation may be a good approximation in 
the IR region. 

We classify all superpotentials W(X) that give rise to confining backgrounds: 
We parametrize their asymptotics for A — > oo as 5 : 

W(X) ~ (logA)* X Q , P,Q e R. (1.11) 

The 't Hooft couplings diverges in the IR as 

p 

g„ / 1 \ 2Q 

\~ E~* Q Hog — j , E^O. (1.12) 

1. Q > 2/3 or Q = 2/3 and P > 1 leads to confinement and a singularity at 
finite r = tq. The scale factor e A vanishes there as 



. (r -r)^ Q > I 

e (r) ~ < r c 



exp 



(ro-r)V(J 3 -!) 



Q 



(1.13) 



2 ' 
3 



where C is a positive constant related to the integration constants. 

2. Q = 2/3, and < P < 1 leads to confinement and a singularity at r = oo 
The scale factor e A vanishes there as 

e A (r) ~ exp[-CV 1/(1 - p) ]. (1.14) 

3. Q = 2/3, P = 1 leads to confinement but the singularity may be at a fi- 
nite or infinite value of r depending on subleading asymptotics of the 
superpotential. 

The above exhaust all cases that confine. All other cases (Q < 2/3, or Q = 2/3 
and P < 0) fail the Wilson loop test. By eq. (|1.7] ), this classification directly 
relates the confining property to the IR behavior of the QCD /^-function. 



4 We assume that the potential, hence the /3-function, do not have singularities at finite A. 
5 Since our results are continuous in the parameters P and Q, our classification also extends to 
any functions W(X) that has a well-defined limit for A — > oo. 
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If Q < 2^/2/3, no ad hoc boundary conditions are needed but UV and IR 
normalizability completely determines the glueball spectrum. This is unlike 
standard AdS/QCD and other approaches. Since the spectrum is completely 
determined from the geometry, and, as discussed above, the latter is in one-to- 
one correspondence with the /3-function, our construction provides a direct 
link between the /3-function and the spectrum. On the other hand, when 
Q > 2V2/3, the spectrum is not well defined without extra boundary con- 
ditions in the IR because both solutions to the mass eigenvalue equation are 
IR normalizable. 6 

For all potentials that confine, the spectrum of ++ and 2 ++ glueballs has a 
mass gap. Moreover, except for the borderline case Q = 2/3, P = 7 , the 
spectrum is also purely discrete. For the + ~ glueballs an extra assumption is 
needed about the strong-coupling asymptotics of the function Z(X) in 

Z(X) ~ \ d , d > 2 as A -f oo. (1.15) 

If this is satisfied, the 0~ + spectrum is also gapped and discrete. We find that 
in QCD d = 4. 

In all the physically interesting confining backgrounds (i.e. those that do not 
require extra boundary conditions in the IR), the magnetic color charges are 
screened. This is shown by studying the Dx branes embedded in the geometry, 
2"U| . This is an improvement with respect to AdS/QCD models, where 



sec 



magnetic quarks are also confined instead of being screened. 

Of all the possible confining asymptotics, there is a unique one that guarantees 
"linear confinement" for all glueballs. It corresponds to the case Q = 2/3, P = 
1/2, i.e. strong-coupling superpotential and /3-function asymptotics: 



W{\) ~ (log A) * A§ , p(\) = -~\ 



} 3 
+ 81ogA + 



;i.i6) 



In this case the 't Hooft coupling diverges with energy as 



3 



3 



A ~ [log^J (1.17) 

in the IR. This choice also seems to be preferred from considerations of the 
meson sector as discussed below. 



6 In an exact solution was studied based on the potential W(X) = Wo + W^Aa in our 
normalizations. This generates logarithmic running in the UV and corresponds to a confining 
background in the IR. However, as it corresponds to Q = 4/3, extra boundary conditions are 
needed at the IR singularity. The spectrum heavily depends on the choice of these boundary 
conditions. 

7 This is also exactly the linear dilaton background in the IR. The spectrum of glueballs is however 
continuous in this case. Because of this, this case is considered no further. 
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The meson sector, assuming Nj -C N c , is implemented by using Nf pairs of 
-D4 — D4 branes embedded in the 5-dimensional background along the lines 



first pointed-out in ten dimensions in [23]. Moreover, the open string tachyon 



field is included and it is dual to the scalar and pseudoscalar quark bilinears 
as first advocated in [F24~fl. In section 5 of 120H we have studied the non-linear 



equation that determines the tachyon profile (vev). It was shown in [23] that 
consistency with the anomaly structure implies that the tachyon field must 
diverge before or at the IR end of space. What we find here is that for the 
confining background the tachyon necessarily diverges at the IR singularity, 
signaling chiral symmetry breaking and the IR recombination of the flavor 
branes. Moreover, it is found that although the tachyon cannot diverge before 
the IR singularity, its derivatives do generically diverge. Such a divergence is 
physically unacceptable. Imposing its absence, determines the quark vacuum 
condensate in terms of the UV quark masses that act as sources for the tachyon 
field. 



As in ||24j| , the spectrum of mesons exhibits (almost) linear confinement due to 
the tachyon potential rather that the graviton-dilaton background. We have 
also studied the masses of the simplest mesons, the vectorial ones which are 
independent of the details of chiral symmetry breaking. They depend in general 
on a combination of the confining QCD scale (appearing in the graviton-dilaton 
data) and the AdS length £. We find that the special background advocated 
above, corresponding to P = 1/2, which provides linear confinement in the 
glueball sector is also the one in which meson masses do not depend on the 
AdS 5 scale, as expected on general grounds. Therefore we obtain a simple 
linear relation between the mass scales of glueballs and generic mesons (this 
does not include the pseudo-Goldstone bosons). 



• We calculate numerically the ++ , 0~ + and 2 ++ glueball spectra for /3-functions 
that interpolate between the standard perturbative QCD regime and the con- 
fining regime, both for the r finite and infinite cases. We compare the glueball 
spectra with lattice results. Although there is no universal consensus on the 
reliability of various lattice results and their relationship to the large- N c limit, 
we attempt a comparison with various confining IR asymptotics. We find that 
the cases of r = 00 are preferred by the data. 



We analyze the axion sector of the theory in section 4 of ]2JJ. We solve the 
equation for the axion to find 



a{r) = {e uv + 2nk) 



dr 



Z(\) 



dr 



(1.18) 



Z(\) 
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where 9jjv is the UV value for the QCD 8- angle and the integer k labels dif- 
ferent large-iV c vacua. In particular we reproduce the well known result that 
the ^-dependent vacuum energy is to leading order in l/N c proportional to 
MirikiOuv + 27r/c) 2 and relate the coefficient (topological susceptibility) to the 
standard QCD /3-function as well as the axion /^-function Z(X). One important 
corollary of our analysis is that a(r), "the effective #-angle" vanishes as a power 
of the energy in the IR, 

0eff(E) ~ £ 3 (log£)3 (1.19) 

This is an interesting result as it indicates that the YM dynamics screens the 
UV 0-angle. 

Before continuing we will comment on what we hope to achieve with this ap- 
proach. As it will become evident, this phenomenological approach does not have 
at this point the status of well controlled approximation. Therefore at best we can 
hope to achieve the following: 

1. To provide a successful and predictive phenomenological model with few param- 
eters. A good analog of such a model is the Lund Monte Carlo that describes 
hadronization in high energy collisions using a basic string model and a dozen 
parameters. 

2. If item 1 turns out to be successful, such a model can provide intuition and 
data, to guide serious searches for constructing a string theory for QCD. 

3. It may provide hints for new phenomena in the theory. A good such example 
may turn out to be our observation that the IR #-angle in large-iV c QCD 
vanishes. 

This introduction summarizes the present paper as well its companion [^(| that 



should be read as a natural continuation of this one. The structure of this paper is 
as follows. 

In the next section we describe in detail the gravitational set up used to explore 
the duals of QCD-like gauge theories. In particular we analyze the general form of 
the scalar potential. 

In section [3], we analyze the equations of motion for the coupled scalar-gravity 
system and derive the precise relation between the full /3-function of the gauge theory 
and the scalar potential. 

In section 4, we derive the UV asymptotics of the solutions close to the AdS 
boundary, using the gauge theory input, namely asymptotic freedom. Then, we 
discuss the qualitative features of a class of IR asymptotics in the deep interior of the 
geometry that lead to confinement. This section also includes a detailed discussion 
of the a 1 corrections near the boundary. In particular, we show that in the expansion 
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near the UV, the a' corrections are always accompanied with the scheme-dependent 
/3-function coefficients. 



Section 5 is devoted to specific examples of geometries that follow from a certain 
choice of the scalar potential. In section 5.1 we present a simple geometry that dis- 
plays asymptotic freedom in the UV and linear confinement in the IR. This example 
is used in [EZ0[] to compute the glueball spectra. In section 5.2 we present a sample 
geometry that displays a Banks-Zaks type fixed point in the IR. 

Section 6 investigates the fluctuations of the various fields in the geometry and 
how the dilaton potential (ie. the perturbative /3-function) modifies the fluctuation 
equations near the UV boundary. 

Various appendices detail our computations. In appendix [A], we present a de- 
tailed analysis of the expansion near the UV. In appendix [B] we analyze a general 
form for the action of a dilaton-axion-gravity system that incorporates the a' cor- 
rections by assuming general forms in the kinetic terms for these fields. We derive 
general solutions of the action and determine the a 1 corrections (in the UV), to the 
various quantities used in the bulk-boundary identification. Appendix |C] describes 
the geometry that follows from the naive effective potential which ignores the a' 
corrections and studies the fluctuations around the AdS vacuum of this potential. In 
appendix [B], we analyze the solutions to the coupled Einstein-scalar system with an 
exponential potential. We derive and classify all the solutions and describe the fixed 
points. Finally, appendix [E] presents another example of a geometry that approaches 
to a conformal fixed point in the IR with an exponential tail in the /3-function. 

As Part II of this work is an essential sequel of this paper we briefly review 
here its structure. It is mainly devoted to the analysis of the non-perturbative 
regime of our construction, i.e. to the IR properties of the 5D geometry. There, 
the general classification of IR asymptotics which lead to confinement is given. The 
qualitative features of the glueball spectra are discussed, and a relation is established 
between the existence of a mass gap and the confining property of the QCD string. 
The dependence of the spectrum on excitation number is discussed, and the IR 
asymptotics that lead to a linear glueball spectrum are identified. The analysis is 
extended to mesons, indicating that the setup provides a concrete realization of the 
holographic implementation of chiral dynamics proposed in J24|. The properties 



of the 5d axion are discussed, and what they may imply for the QCD ^-parameter. 
Finally, numerical computations of glueball spectra in concrete models are performed. 
The models are defined in terms of an exact /3-function that interpolates between the 
desired UV and IR asymptotics. The detailed structure of the Part II can be found 
in its introduction. 
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2. Motivating the gravitational dual of a large-N c gauge the- 
ory 

The gauge theories of interest in this paper, are four-dimensional U(N C ) gauge the- 
ories at large N c . In particular, we assume the presence of no further adjoint fields, 
and therefore, the holographic dual theory is expected to live in five dimensions. 
Fundamental matter can be present but we will assume here that the number of 
flavors Nf <C iV c . Therefore, quarks can be eventually incorporated as four-brane 
probes inside the five-dimensional geometry. 



2.1 The spectrum 

The relevant non-critical string is therefore five- dimensional 8 . As there are no fermionic 
gauge-invariant operators in pure YM theory we do not expect the associated string 
theory to have space-time fermions, (No NSR and RNS sectors). 

^From the NSNS sector we will obtain as low lying fields a metric g^ u , a dilaton, 
and a two- index antisymmetric tensor B^. The generic bulk fields should be in one- 
to-one correspondence with the low-dimension gauge-invariant operators of the gauge 
theory. In particular the (classically) traceless stress tensor T^ u = Tr[F 2 u — \ri^ u F 2 } 
is dual to the five- dimensional graviton g^ u . Tr[F 2 ] is dual to the dilaton 0. B^ u is 
expected to be dual, in analogy with M = 4 sYM, to a higher dimension operator 



of the form Tr[d a b c F a F F c ] [3d]. As it will be trivial in the vacuum solution, due 
to Lorentz invariance, we will neglect it in the sequel. It is however expected to 
generate a tower of 1 h glueballs in the theory. 

The RR sector massless fields are summarized into the tensor product of two 
5-dimensional spinors. The product of two five- dimensional spinors gives the fields 
strengths F , F 1 , F 2 , F 3 , F^, F 5 . There is an automatic duality condition for this prod- 
uct that relates F 5 ~ F , F 4 ~ F 1 and F 3 ~ F 2 . Therefore only F ,Fi,F 2 are 
independent. 

F\ = dCo generates an axion field Cq = a that is dual to Tr[F A F] in the YM 
theory. The dual of the axion field-strength is a four-form field strength F^ and the 
associated three- form C 3 couples to domain walls (that separate different fc- vacua). 

F 5 = dCi ~ F generates a four-form that seeds the D 3 branes responsible for 
the U(N C ) gauge group. Its dual is a zero-form field strength that couples to bulk 
instantons. 

Finally F 2 = dC\ generates a vector and its presence seems puzzling as there is 
no candidate YM operator that it is dual to. It turns out however that such a form 



couples minimally to the baryon density on flavor branes | 2"3[] , and its presence is 
therefore justified. It is certainly trivial in the QCD vacuum as it will break four- 
dimensional Lorentz invariance otherwise. Its kinetic term is suppressed by as 



For a similar discussion of the spectrum, see 25 
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that of the axion. It can however be used to turn non-trivial baryon number densities 
in QCD. It is expected to generate a tower of spin-1 glueballs comparable in profile 



to the + ones. We will analyze it further in p0 |. 

We would like at this point to stress that the presence of RR fields is important 
with matching with the behavior expected from QCD. In particular, the special 
dependence on the dilaton of terms in the effective action containing RR fields (in 
particular the axion), is in complete agreement with the large N c suppression we 
expect for such terms in the gauge theory. 

Flavor is expected to be generated by space filling D 4 + Z) 4 brane pairs filling 
(most of) the five dimensional bulk. They couple individually to a C5 form. Overall 
neutrality is needed however as it equivalent to anomaly cancellation in the YM 
theory. 9 As the branes are space-filling, the coupling to C5 will not play a role. 
Subleading couplings to other RR fields are essential though and will be discussed 
in 



The structure of the RR sector is similar to what we would obtain from an 



orientifold of the OB theory [42. 43. 16] compactified to 5 dimensions. Several authors 
advocated the relevance of a closed string tachyon in the vacuum structure of QCD 
motivated by the ten dimensional type-0 string theory. When Nf = we see no place 
for such a closed string tachyon, as there is no such low-dimension gauge-invariant 
operator in the gauge theory. If type-0 theory has any connection to pure YM 
theory, then the ten-dimensional tachyon must become rather massive in the five- 
dimensional theory. Of course if we add Nf flavor branes with Nf <C iV c then the 



open string tachyon is essential for understanding meson physics, |24| . However, the 



tachyon profile does not back-react to correct the gauge theory vacuum in this case. 
A (closed) string tachyon may be relevant in understanding the vacuum structure 
of QCD in the Veneziano limit, Nf ~ N c . In this case it is the avatar of the open 
string tachyon of the flavor branes, which, because of their large number is becoming 
effectively a closed string state (in the sense that it affects the vacuum structure at 
leading order in 1/N C ). 

To summarize, when Nf = the relevant propagating bulk fields are 5^,0, a, 
while there is also a non-propagating F5 field strength. When N c ^> Nf 7^ 0, there is 
also an open string tachyon T (together with gauge fields on the flavor branes), which 
is a Nf x Nf complex matrix charged under gauge fields of the U(Nf)i, x U(Nf)n 
symmetry, [24]]. In this paper we will not consider the dynamics of flavor. But we 



will include the contribution of the tension of flavor branes in the effective action 
in this section only to indicate the scaling of various terms in the effective dilaton 
potential. 



9 This appears to be rather general in holographic setups. Cancellation of bulk RR charges is 
equivalent to the absence of gauge anomalies of the large- N c boundary theory. 
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2.2 The associated branes 

There are several potential branes associated to the various bulk form fields that exist 
in the theory. The RR four-form couples to D 3 branes, responsible for generating 
the gauge group. 

The RR zero form, dual to the YM #-angle couples electrically to D_i instantons, 
which should be indeed the YM instantons. It is interesting to mention here that 
such bulk instanton solutions in the AdS^ case were identified with the standard 
YM instanton solutions with the dependence on the radial holographic coordinate 
playing the role of instanton size. In our case we expect that small instantons, as 
they will be close to the AdSs boundary, resemble almost accurately the standard 
YM Instantons. Large instantons however, feel the non-trivial IR geometry and 
their form is expected to part substantially from the standard YM form. It may be 
interesting to study them as they might indicate the correct instanton measure in the 
IR of YM theory. This should be the five-dimensional volume form in the non-trivial 
metric that describes the YM vacuum, together with a power of the dilaton that 
needs to be determined. 

The magnetic duals of the _D_i instantons are D2 membranes and they couple 
minimally to the dual C3 form. They are domain walls in four-dimensions separating 
different fc-vacua. A point-like D 2 generates a monodromy for the bulk axion field 
in the transverse two dimensional space. 

There are also point-like D branes coupled to the RR one-form C\. They are the 
baryon vertices which can tie up N c quarks to form a large- N c baryon. Their duals 
are D\ strings that couple to the dual two-form C2, which represent the flux tubes 



stretched among magnetic quark sources. They will be used in |20| to investigate the 
interaction between magnetic quarks in confining vacua. 

Finally in the NSNS sector we have NSo branes that are five- dimensional duals 
of the fundamentals strings. They are charged magnetically under the NSNS two- 
form B^, and couple directly to the flavor U(l) gauge-fields. The S-duality map 
would suggest that they can be viewed as vertices for magnetic baryons. This is in 
agreement with the fact that such magnetic baryons are expected to have masses 
that scale as 0(N%) compared to ordinary baryons whose masses are of order 0(N C ). 
The reason is that they contain N c magnetic quarks, each having a mass of order 



10 



2.3 The effective action 

The (minimal) two- derivative effective action of the perturbative string theory is 



S 5 = M» J d^Vg [.-* (r + 4(fl«» + |) - ^Fl - if? - 2f 



e 



(2.1) 



10 See however a slightly different interpretation in [E5j 
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where as usual F\ = d^a and we have used the 5-form instead of the zero form. The 
last term is due to space filling D4 — D4 brane pairs, with Nf is proportional to the 
number of flavor branes. We have explicitly indicated the string scale £ s and the 
Planck scale 



1 



5c = 10 - D 



Anticipating the connection with the gauge theory let us define, 



(2.2) 



A = N r e 4 



(2.3) 



4 p 

Passing to the Einstein frame by g^ u = \^gj^ v we obtain: 



M 



d 5 x^/g 



'^A) 2 5r . 



(2, 



We are now at the point where we can solve for the five-form (that controls the 
number of D% branes). Starting from the Einstein frame action 



6S* 



M 3 
2^5! 



dPxy/g A-t F 2 , 



(2.5) 



the equations of motion are d*(X 3 F 5 ) = with the solution 



N r , 10 „ iV r , 10 e,,, ...... 



(2.6) 



where E is the totally antisymmetric tensor and N c is dimensionless and propor- 
tional to the number of colors. This is the definition of N c , and it is obviously 
unormalized, unless we know explicitly the D$ brane solution of the associated string 
theory. Inserting back in the equations of motion and from this reconstructing the 
action, we obtain 



J d 5 xy/g 



e 3 ' 



so that 



with 



S 5 = M 3 N 2 J d 5 xy/g- 



R _Hd\y 



A 2 



3 A 2 



27V 2 



(da) 2 + V{\) 



V(X) 



At 
I 2 



1 



5c — x\ A 2 

2 



Here we defined the ratio of number of flavors over the number of colors, 



(2.7) 



(2.1 



(2.9) 



x = 



(2.10) 
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This potential is analyzed in detail in appendix |C[ It has a minimum at a finite 0(1) 
value of the t 'Hooft coupling 11 1273 . 



This gives rise to an AdS 5 solution that was found and analyzed in ||, where 
it was conjectured to describe the Banks-Zaks fixed points (for Nf ~ N c , i.e. x ~ 
1). One might think that a dilaton flow around the AdS^ solution might describe 
the gauge theory physics we are seeking. A study of the dilaton perturbations in 
appendix shows however that when Nf <C N c , the dimension of the operator dual 
to the dilaton perturbation is well above 6, and this does not compare well with the 
dimension-4 perturbation we are looking for. On the other hand, it is expected that 
for Nf ~ N c , the solution will be importantly affected by the tachyon field. Therefore 
an analysis in this case that neglects the tachyon effects may be unreliable. Thus 
it is plausible that the AdS 5 solution found is relevant for gauge theories in their 
Veneziano limit, as argued in [§]. 

Another problem of this AdS vacuum is that, as we show in appendix |Cj, the t 
'Hooft coupling at that point is of order one, and this is not something we expect in 
the UV region of QCD. 

We are therefore led to the suggestion that the UV of QCD should be described by 



the — > — oo asymptotics of the potential in ( |2.8|) . To test this claim, we must study 
the asymptotic solution of the classical equations and compare it to expectations 
from QCD. To do this we make a Poincare-invariant ansatz for the metric. The 
equations of motion in the covariant form are, 



F _t 

u 3 



with the Einstein tensor defined as, 



= R^ u — -^Rdfiu, (2-12) 

We neglect the QCD axion a in our system as it is not expected to affect this 
solution. 12 

Let us digress here to introduce the following two coordinate systems we use in 
the paper. The "conformal" coordinate system is given by, 

ds 2 = e 2A{r) [dr 2 + ru j dx i dx 3 '] , (2.13) 



n We stress again that this is the definition we use for the 't Hooft coupling and it is unormalizcd. 

12 In large- N c theories with strong IR dynamics like QCD, Witten has argued long time ago p6| , 
that although instanton effects seem negligible at large--/V c , there is non-trivial 9 dependence in 
QCD dynamics. This is in particular responsible for the resolution of the rj mass problem via the 
Witten- Veneziano mechanism, f27f . Although we expect via such considerations a to have a non- 
trivial profile in five dimensions, its "order of magnitude" is 0(1/N C ) and therefore its back-reaction 
to the other bulk fields can be neglected to leading order p8[. 
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and the "domain wall" coordinate system is given by, 



ds 2 = e 2A{u) [r^cfeW] + du 2 . (2.14) 

The two are related by 13 , 

du = e A{r) dr. (2.15) 

The advantage of the conformal variable r is that in an asymptotically AdS back- 
ground, it is directly related to the energy of the gauge theory as E = 1/r in the far 
UV. Also the corrections to the AdS metric that we describe in this paper take the 
simple form of logr. On the other hand, the domain-wall coordinate u facilitates the 
solution to the differential equations and uncovers the relation between the integra- 
tion constants and the related gauge theory quantities. Henceforth, we will denote 
by prime, a derivative with respect to u and by dot, a derivative with respect to r. 
The equations of motion in the conformal coordinate read 

12v4 2 — —(f) 2 — e 2A V = , QA + QA 2 + —cj) 2 — e 2A V = (2.16) 
3 3 

+ 3i0+^| = O (2.17) 

We shall also need the Einstein's equations and the scalar field equation in the 
domain- wall coordinate, in the following form: 

A n = W_ + y_ £ii _ _!(0')2 (2.18) 
9 12 ' 9^ 1 ' V ' 

(2 - l9) 

What type of solution do we expect in the UV regime (r — > 0) (or u — > — oo)? 
As QCD becomes scale invariant at E — oo, we expect the space to asymptote to 
AdS*,. Moreover as j- ~ log-E we expect that ~ — logr. Therefore, the geometry 
is expected to be AdSs modulo logarithmic corrections. This asymptotic expansion 
is described in detail in appendix 0. 

The potential in ( |2.8| ) as — ► — oo, behaves as V ~ e a< ^ with a = |. In appendix 
we analyze in detail the solutions to the single exponential potential. We show 
that in such cases the asymptotics of the metric is quite different from the expected 
behavior outlined above. Indeed one finds, as r — » 0, 



e 



16 . 

A(r) _ J r^ra if a ^ | 



e r if a = |. 



13 For simplicity we denote the warp factor in both of the coordinate systems by A i.e. A(u) 
A(r(u)). 
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In fact, the dilaton potential that has the desired behavior in the UV, must 
asymptote to a (positive) constant, as — > — oo. The natural question is whether 
such a constant can appear from string theory considerations. 

In five dimensions there are higher-order a' corrections, that contribute to the 
dilaton potential. This may sound implausible, as such corrections carry a number 
of derivatives and are not therefore part of the potential. However, it should be 
remembered that the bulk theory contains the (non-propagating) four- form C4, which 
is responsible for the charge of D% branes sourcing the bulk solution. Its field strength 
is F5. At the quadratic level, it generates the term in the potential, upon 
dualization. It goes without saying that its higher-derivative avatars in the effective 
action will contribute several more terms in the dilaton potential. For example the 
n-th order term (on the sphere) in the a' expansion is 

m 3 r 

S "~p^J d 5 xe 2 ^^ (F 5 ) 2n (2.20) 

where the dilaton dependence is dictated by the special transformation properties 
of RR field strengths. As explained in detail in appendix [B], upon dualization, such 
terms translate into a dilaton potential of the form 

±0 OO 

V F5 ~ e —J2 a n(N c e*r n (2.21) 

where the a n are dimensionless numbers reflecting the coefficients of the higher deriva- 
tive terms like Q2.20| ). Note that we included the central charge deficit as the first 
term in this potential. 

This potential is remarkable for two reasons. Firstly, it is a function of the 't 
Hooft coupling A ~ N c e^, and remains finite in large-iV c perturbation theory. The 
second is that it has a regular expansion for small 't Hooft coupling. On the other 
hand it does not contain the term we are looking for, namely a (positive) constant. 



In appendix p.3| we indicate how such an effective term, generating an asymptotic 
AdS regime can be generated by higher curvature corrections. We will not however 
pursue this line of thought here, instead we will add such a term to the dilaton 
potential by fiat, as was also done in [|X2j] for example. 

The next question to address is what the other terms in the potential should be 
in order to reproduce the expected logarithmic running of the gauge theory coupling 
constant A, that is compatible with QCD perturbation theory. To answer this ques- 
tion, one must identify the UV energy variable, and this is essentially the coordinate 
r up to logarithmic corrections. Matching then the expected equation 

<IX (3{X) = -b X 2 + M 3 + M 4 + ■ • • (2.22) 



dlogr 
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with the bulk equations of motion in (|2.16|) , ( 12.17 ) we obtain that the weak-coupling 
expansion of the potential should be 



V(\) = ^T / V n \ n (2.23) 



n=0 



To leading order, A ~ N c e^ as can be ascertained by the world-volume couplings of 
D3 branes 14 . Higher a' corrections involving the five- form on D3 branes can change 
this identification to 

00 

X = J2 c n (A^) (2n+1) (2.24) 

n=0 



as is shown in appendix B.2. In case of M = 4 YM, the leading coefficient is 



Co = 47r. However, in general we consider a possible multiplicative renormalization in 
the identification of A, hence we keep Co arbitrary. Moreover, changing Co is equivalent 
to changing 6 the first /3-function coefficient. Therefore we will eventually set Cq — 1 
and we will allow ourselves to vary 6 to reflect this normalization ambiguity. 

Such corrections affect the /3-function beyond the first two non-trivial orders, 
(associated to bo, b\). Since higher coefficients are non-universal, we choose to neglect 
the renormalizations in fl2.24| ), and set c n = for n > 1. In the analysis of appendix 



[7a| we have set Cq to 1. Reinstating it here we obtain 



a 23 - 36& 324^ + 124 + 189^ 

1 + g(M) + 34 ° (M) 2 - 2^ — - 7 ^(M) 3 + 0(X 4 



12 



a 23 - 36^ 324^ + 124 + 189^ 
1 + g(&oc iV c e^ + -^A(6 oCo iV c e^) 2 - 2 * b <L(b,c Q N c e*Y + 

(2.25) 

where I is the AdS radius. We conclude that without knowledge of cq, the effective 
one-loop term b e ^ = b c is unknown. 

In appendix it is shown that changing the constant 60 amounts to a redefinition 
of the perturbative QCD scale A as well as a redefinition of the non-universal coeffi- 
cients of the /3-function. This is nicely reflected in the fact that the dilaton potential 
( |2.25| ) depends only on the invariant ratios bi/b Q +1 . Moreover as shown in the same 
appendix, the first logarithmic correction to the UV AdS metric is universal, and in 
particular insensitive to either bo, or cq. In view of the above, we will postulate from 



now on, a dilaton potential with a weak coupling expansion as in (|2.25 ) with cq = 1. 



14 If the space has more than five dimensions or there are other bulk scalar field involved, this 
identification can change. However this is not the case here. Another issue involves finite renor- 
malizations of the 't Hooft coupling identification. Such corrections cannot be excluded but they 
cannot affect the qualitative analysis, as they involve a constant rescaling of the dilaton. We choose 
therefore to ignore them. 
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3. Holographic large N c YM: the vacuum structure 



Before we embark on a study of the solution to the equations of motion that describes 
the vacuum structure of the theory, we review the action we are going to use 



S = M 3 N 2 J d 5 x^/g 



4(«9A) 2 



Z{\) 
2AT? 



(da) 2 + V(X) 



(3.1) 



with 



V(\) 



8 236g - 36b! a 
9 M + 81 A + 



(3.2) 



and where we have added a general kinetic function Z(X) in front of the axion. As 
the axion kinetic energy is suppressed by 1/AT? it will be neglected in the following 
investigation of the vacuum solution. In the action (|3.1| ), N c appears only as an 
overall coefficient, i.e. it only changes the definition of the 5-dimensional Planck 
scale. Therefore in the large N c limit, quantum corrections can be safely neglected, 
as in the standard, critical AdS/CFT correspondence. All kinetic and potential 
terms are independent of N c when written in terms of A. Therefore, it is useful to 
define an "effective" dilaton field by 



$ = log A 



+ logiV c 



(3.3) 



From now on, this is the field we will use, and we will refer to it as the "dilaton" 
without any further specification. 

3.1 General geometry and the geometric invariants 

In this section we introduce a rewriting of the equations of motion that facilitates 
the construction of the desired 5D geometry and its connection to the gauge theory. 
We allow for a general dilaton potential V(<&) in (|3.1| ) that admits an expansion of 



the form (^2) for small A. 

We start from the domain- wall parametrization of the metric in (|2.14j) . We use 
the Einstein's equations fl2.18|) and the scalar field equation ( |2.19| ). It is a well- 
known fact that these second-order equations can be reduced to a system of first- 
order equations 15 . This is generally done by introducing a superpotential 16 . Here 
we find it more convenient to rewrite the equations in terms of the derivative of the 
superpotential with respect to the dilaton. For this purpose, we introduce the phase 
space variable X: 

(3.4) 



X 



3A' 



15 This well-known technique in differential equations was discussed in this context in |29|| . 
16 It goes without saying that this amounts to choosing one of the two initial conditions for the 
dilaton. 
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This variable is useful in our set-up as it is closely related to the gauge theory (3- 
function. 

In terms of X, the second order system can be rewritten as the following first- 
order system of equations, 

($'f = !y xVl^ x# (3.5) 

(A0 2 = §e-i/-*-*» (3.6) 
dX , d\ogV,(X 2 -1) . „. 

Here, V > is the asymptotic value of the dilaton potential as $ — > — oo. Given 
X as a function of one solves (|3.5|) to determine the dilaton and (|3.6|) to obtain 
the geometry. Our strategy in solving the system of equations is essentially reversing 
( |3.7|) and determining the potential V(<&), given the function X($). In the next 
subsection, we explain how to determine X, by using the gauge theory input. It is 
completely determined by the full /3-function of the gauge theory. We analyze certain 
fixed points of ( p.7| ) in the appendix |D[ 

Here, we assume that X is known, and calculate the invariants of the geometry. 
This is necessary in order to understand the singularity structure of the solution in 
the IR. 

We solve ( |3.7|) for V as follows: 

V($) = V (l-X 2 )e-*f-°° Xd *, (3.8) 

where Vo > 0. 

The leading scalar invariants, that are constructed from the Einstein tensor (de- 
fined in Q2.12Q ) are, 

E = E/ = 30(A') 2 + 12A", (3.9) 
E^E^ = 36 (5(A') 4 + 4A"(A') 2 + (A") 2 ) . (3.10) 

Using the equation of motion (|2.18|) together with (|3~5| ) and (|3.6| ), we obtain the 
following results for the first-order invariant: 

E = V (^-AX 2 )e-l^ Xd ^ (3.11) 

We show in the next section that the presence of an asymptotic AdS geometry near 
the boundary requires X — > there. From ( |3.11| ) we observe that E — > 5Vo/2 
there, as expected from the AdS geometry. This provides a consistency check of our 
formulae. The second order invariant is, 

E-E=^-(2- 11X 2 (1 - X 2 ))e~f f-~ xd *. (3.12) 
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We also present the relation of the phase variable and the superpotential. In 5D, 
the superpotential is obtained from the scalar potential by solving, 

3\ 2 fdW\ 2 /3 x3 



w -\i) [w) =U> (3 ' 13) 

The first order equations that follow from this superpotential are, 

A = - A -W , *' = ^. (3.14) 

9 ' d$ v ; 



Comparison of (|3.14|) with (|3.6| ) and (|3.5|) shows that the superpotential is given as, 



W=^e-i£- x » (3.15) 
8v 3 

It is also useful to note that X is proportional to the derivative of the logarithm of 
the superpotential: 

4 rf$ v ' 

3.2 The phase variable X and the /5-function 

In this section, we show that given i) the relation between the 't Hooft coupling A and 
4> and ii) the relation between energy scale of the gauge theory, and the radial variable 
of the 5D geometry, one can relate the phase space variable that we introduced in 
(|3.4D to the /3-function of the gauge theory. Both relations receive a' corrections, but 
we shall first ignore them and consider their effects in section |4.2| . 

To fix the first relation, we consider the metric ( 2.14Q in the Einstein frame. 



There are at least two arguments that lead to the energy-radius dependence: The 
first is the Polchinski-Strassler type of argument that uses the gravitational 
red-shift to relate 

E = v^oo E bu ik, (3-17) 

where the LHS is the energy in 4D field theory measured by an observer on the 
boundary of the 5D geometry and the RHS is the energy of a point-like excitation 
of the gravitational system. This gives, 

log(.E) oc A. (3.18) 

The second argument is to make an infinitesimal shift u —>■ u + 5 in the metric and 
to observe that this shift can be produced by a dilatation x M — > (1 + p)x tl . This 
dilatation changes the 4D energy by, log(-E) — > \og(E) — p. Equivalently, one makes 
the above shifts in u and x and demands that the metric does not change to first 
order in p. This fixes, 

p = -A'5 
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hence one arrives at the same result (|3.18 



We note that the relation ( 3.18|) is the most natural in an asymptotically AdS 



geometry where A — > —u/l as u — > — oo and one identifies logE = u in the original 
AdS/CFT correspondence (|£§). However, we stress that any function whose leading 
behavior agrees with ( |3.18| ) in the UV is a candidate definition for the field theory 
energy. 

Another necessary condition is to have a monotonically decreasing function of 
the radial variable u. This is necessary in order to have a well-defined direction for 



the renormalization group flow. From the Einstein's equation ( 2.18|) , it is clear that 



( |3.18|) is guaranteed to satisfy this criterion. This is because, it follows from Q2.18j ) 



that A' is a monotonically decreasing function of u. Since it is already negative at 
u = — oo, A' is forced to be negative in the entire range of u. Thus, A itself should 
be a monotonically decreasing function of u. 

We should also note that the scale factor of the metric in the Einstein frame is 
expected to be associated to the C-function of the associated theory. In particular 
with our definition we expect that E 4 ~ e 4A is proportional to the number of degrees 
of freedom of the theory with energy up to E. 

Now, we come to the identification of the /^-function. As discussed in Section 0, 
we identify the 't Hooft coupling as: 

A = N c = e*. (3.19) 

where we chose to set the coefficient cq to one here by reabsorbing it in the potential. 
On the other hand, the /3-function of the gauge theory is given by, 

i3(A > = SktfEj = TA < 3 ' 20) 



The definition of the phase variable X in ( p.4|) together with the gauge theory- 
geometry relations ( |3.18j ) and ( |3.19| ) leads to the following identification of the bulk 



and boundary quantities: 

X = m. (3.21) 
3A v ' 

In the next section, we shall make use of (|3.21|) in order to formulate restrictions on 

the geometry that follow from the desired properties of the gauge theory. 

4. Asymptotics of the geometry 

4.1 Asymptotic freedom and the UV geometry 

In an asymptotically free gauge theory, the 't Hooft coupling behaves as, 

A oc - — , for large E. (4.1) 
logE 
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Using our identification of the energy scale of the gauge theory with the scale factor 
of the metric in the Einstein frame, i.e. logE oc A, we learn that this requirement 
means, $ oc — log A, or in terms of the phase-space function X, 

X oc -- oc -A. (4.2) 
A 



Making use of (O) in (|3.8|), one learns that the geometry dual to an asymptotically 



free gauge theory, requires the following form of the scalar potential near A = 0: 

V — V + ViX H (4.3) 

as was already announced in section |^. Therefore asymptotic freedom in the UV 
requires an asymptotically AdSs geometry. In addition to this, it requires an expan- 
sion of the form (|4.3| ) around A = 0. In particular, the presence of the first order 
term V\ is crucial for the (leading) logarithmic running of the coupling constant. 

One can easily make a consistency check on the condition, (f4.2|) . Having an AdS 
geometry in the UV gives A oc u near the boundary (u — > oo). Therefore, near the 
boundary we have, 

A oc -. (4.4) 

u 

Using the eq. ( |3.6| ) or ( |3~5| ) we see that this results in the following asymptotic form 
for X: 

X oc A + C(A 2 ), as A 0. (4.5) 



In particular, as discussed in section ||, the simple set up in (|2.4| ) does not satisfy 



this asymptotics. Although we presented ideas on how the two requirements can be 
reconciled, in this work we will assume the presence of the suitable dilaton potential. 

4.2 The a' corrections 

We now generalize the above discussion by taking into account the possible a' cor- 
rections near the boundary. As is clear from ( |3.8| ) the requirement ( |4.3| ) follows from 
requiring Q4.5| ). We would like to determine the asymptotics of X by taking into 
account the possible corrections that fall into the following three classes. First, there 
is a subclass of the derivative corrections to the Einstein-Hilbert action studied in 
appendix ([B]). Second, there are the derivative corrections to the probe D 3 brane 
action studied in appendix (|B.2|) . The combined effect of these corrections can be 



stated as, 

A = A, (l + Cl X 2 s + (2a 2Cl + | + cl)Xt + 0(X 6 S )) . (4.6) 

Here, 

A s = e*N c , (4.7) 

is a purely bulk quantity, while we use A to denote the 't Hooft coupling of the gauge 
theory. The unknown coefficients cij and q are defined in ( |B.2| ) and ( B.40| ). 
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The third type of correction involves the definition of the gauge theory energy 
scale in terms of the metric scale factor. Using the fact that all of the functions in 
the geometry can be expressed in the variable A s , and that, in the far UV, the radial 
variable in AdS exactly corresponds to the energy scale, we may write, 

^ = /(A.) = 1 + fiXl + f 2 Xi + ■ ■ • (4.8) 

with fi > 0, for all i. This expansion anticipates that the variable X s asymptotes to 
zero in the UV, and this will be checked below. We note that a specific choice for 
the function /(A s ) corresponds to a specific scheme 17 in computing the /3-function of 
the field theory. The reason that the expansion variable is A 2 instead of X s is due to 



the form of the a'-expansion, as derived in Appendix [B.3 
Now, we can connect X to (3: 

P = d^E = ^ S JX^E = Z\7 3XAJ(As)_1 - (49) 
Eqs. (p|) and (g~JD give, 

X = ——A + —A 2 + (— - A - Cl b + ^ V 3 + f- + - —) # (4-10) 

3 3 \3 3 3 J \3 3 3 J K ' 

+^ (6 4 + &2(4ci - /1) + 2ci/i6o + ci&o - /1&0 + /260 - 4a 2Cl 6 - 2c 2 6 ) A 5 H 

We remark that the first two terms are completely determined by the scheme indepen- 
dent gauge theory input, i.e. the first two /3-function coefficients, 60 and &1 18 . A very 
interesting corollary of (|4.10|) is that the unknown coefficients of the a' expansion, 
namely a i: Ci and fi are always accompanied by the scheme dependent /3-function 

coefficients b 2 , 63, This, for example, makes it possible to set all the fi to zero by 

choosing an appropriate scheme in the computation of the full /^-function. In what 
follows we shall assume that this scheme is chosen and we shall set fi = for all 
i. We call it the "holographic scheme". A final observation in Q4.10Q is that the a' 
corrections to the effective action discussed in Appendix [FJ, (associated to a,), enter 
first at the fourth order in the expansion. 

So far we have dealt with the a' corrections near the boundary. The possible 
corrections in the IR regime shall be discussed in section 

4.3 The UV geometry 



We can obtain the UV asymptotics of the metric by using ( |4.1(J| ) and the formalism 



of the previous subsection. The asymptotic form of the metric to leading, subleading 

17 The scheme dependence of the renormalization group running in the holographic setup has been 
discussed in 

18 As explained earlier, in the correspondence of the bulk and the boundary quantities, neither 60 
nor bi but only the ratio b\jb\ has an invariant meaning as we do not know the correct normalization 
of A 
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and next to subleading order, turns out to be independent of the corrections discussed 
above: 

ds 2 = e 2 Tul +2b ^ u+ -dx 2 + du 2 , (4.11) 



where I is the radius of AdSs: 

£2 _12 

v - 

We have defined the following ratio of the /3-function coefficients, 

&i 



(4.12) 



b = -± (4.13) 

This ratio has the interesting property that it is left invariant under a multiplicative 
renormalization of the coupling constant. This follows from the Gell-Mann-Low 
equation and explained in detail in appendix [A[ 

On the other hand, the asymptotic form of the dilaton to subleading order is, 

b X = - + (b-^)logu + --- (4.14) 
u 9 

One can also easily obtain the asymptotic form of the Ricci scalar and find that 
R — > —20/£ 2 : We have the AdS§ geometry in the UV and the subleading term in 
( ft.llp can be viewed as a perturbation on AdS 5 . 



We also present the same asymptotics in the conformal r-coordinate system: 



2 V 91og(rA) 9 log(rA) 2 

A _ 1 , b log(-log(rA)) 

6 log(rA) + 6 (log(rA)) 2 + "' 1 ' ] 

where I is the AdS radius and A is an integration constant, that corresponds to the 
perturbative QCD scale in the dual theory. 



We remark that the leading and the subleading terms in the metric, ( 4.11 ), de- 
pend neither on the energy-radius coefficients fa, nor on the particular /5-function 
coefficients 6j. Therefore the first terms in (|4.11|) are universal and solely follow from 
the requirement of asymptotic freedom. However, the subleading term of the metric 
already introduces the invariant QCD scale A. Of course the information on the 
/^-function is encoded in the leading term of the dilaton solution. 

4.4 Confinement and the IR geometry 

Given the geometry, a natural question to ask is whether or not this geometry de- 
scribes quark confinement in the IR. On the gauge theory side, a criterion for con- 
finement is that the Wilson loop operator exhibits an area law behavior at large 
distances. This corresponds to an asymptotic linear potential between two external 
quarks. In the 5d dual this criterion can be tested by calculating the classical action 
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of a fundamental string world-sheet that extends to the Wilson loop at the AdSs 
boundary, gued in |34], [35]]. 



A detailed analysis with a classification of all the confining backgrounds is pre- 
sented in Here we would like to discuss the qualitative features of the set-up by 
focusing on one class of asymptotics: 

1 



X(X) 



log A 



A -> oo (4.17) 



The asymptotics of the geometry can be obtained using the definition ( |3.4|) : 

A - A ^iL (418) 



and $ from (|3.5| ). In the conformal variable, one obtains the following IR asymptotics 
(as r — > oo) 

A = -Cr a + subleading, (4.19) 
where C is an integration constant and we have defined, 



The dilaton behaves as, 



a = ^Ta (4 ' 20) 
3 

$ = -Cr a + subleading (4.21) 

It is shown in [^tj that these backgrounds confine for a > 1. To see this, one needs 
to take into account the subleading term in ( |4.21|) . In the borderline case, a = 1, the 
geometry is interesting and well-known. In the string frame the metric asymptotically 
becomes, 

ds 2 -> dx 2 + -d$ 2 , $ -> oo. (4.22) 
9 

This is the linear-dilaton background. The Ricci scalar in the string frame therefore 
vanishes in the IR. In the Einstein frame, the metric has a singularity at a finite 
value u = Uq (corresponding to r — > oo. There, the scale factor exp[A] shrinks to 
zero. This is due to the fact that e* and therefore the gauge coupling diverge at that 
point. The stringy corrections to the geometry in this regime can be neglected as in 
the string frame the background is known to be exact to all orders in a'. 

The linear dilaton asymptotics are the borderline of confining backgrounds. The 
spectrum of mesons and glueballs is obtained by considering the fluctuations of the 
geometry. Although exhibiting a mass gap, the spectrum is partly continuous for 
a — 1. Therefore, we discard the borderline case a = 1 (a = 0) for physical reasons. 
In the case a > 1 however, the spectrum turns out to be discrete. 

For a > 1 (or a > 0) the linear dilaton background in (|4.22j) is modified into the 
following asymptotic metric: 

ds 2 -> $^dx 2 + ^$-Trd$ 2 , $ -> oo. (4.23) 
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The Ricci scalar in the string frame reads asymptotically , 

R s Aa(7a - 3)$t a - 2 . (4.24) 

We observe that for the confining backgrounds a > 1 (or < a < 3/4) the string 
frame Ricci scalar vanishes in the IR. This implies that the singularity in the string 
frame is only due to the diverging dilaton. This is a good type of singularity where 
the string moving in this background feels very small curvature corrections. 

The effect of a' corrections to the confinement criteria is an important question 
to which we turn. Fortunately, as happens with "good" holographic singularities, 
the low-lying fluctuations have an effective potential that shields them from the 



singularity. Therefore, the a' corrections are under control in this case, see [20 . 

Let us finally note the behavior of the dilaton potential and the various geometric 
invariants in this region. The dilaton potential follows from 



V « ^oAt(logA)^ (l + O(j^)) • (4-25) 

Similarly one obtains the following asymptotics for the geometric invariants in the 
Einstein frame, 

E ~ R~ (<9$) 2 ~ A3 (log A) ^ (4.26) 

We see that this is the linear dilaton background with logarithmic corrections. 
In the a = 1 case the leading behavior of the geometry is exactly that of the linear 
dilaton background. 

5. Examples of background geometries 

In this section we present examples of full background geometries that exhibit the 
desired UV and IR asymptotics discussed in the previous two sections. The first is 
an example showing asymptotic freedom in the UV and confinement in the IR. For 
comparison, we present an example which is asymptotically free in the UV, but has 
a conformal fixed point in the IR. Finally, in appendix [E| we present a third example 
which is asymptotically free in the UV and has a /5-function that asymptotes to zero 
in the IR with an exponential tail. 

5.1 Standard QCD type two-loop /3-function 

We start from the "exact" /3-function: 

«A) = -"5*1 T - (5.D 



3 + 26 A 



! o A 3(1 + A 2 ) 1 + (2ft °+f log(l + A 2 

with &o,i > 0, a > 0. This expression is chosen so that one has the desired asymptotics 
in the UV and IR and such that there are no poles or branch cuts in A. In particular at 
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V 




Figure 1: Dilaton potential as a function of A for the solution studied in section 5.1. 



weak coupling it reproduces the QCD result to two loops, (3 = —b X 2 — bi\ 3 + C(A 4 ). 
At strong coupling it behaves as (3 = — |A + 3a + ■ ■ ■ as required by ( |3.21| ) and 
( f4.17| ). It also avoids unwanted singularities or branch points at finite values of A. 



According to ( |3.8[ ) the dilaton potential that corresponds to this specific (3- 
function is, 

V(A) = Vb(l - A(A) 2 )(3 + 2b \f* (l8a + (26 2 + 3&i) log(l + A 2 )) ^ , (5.2) 

where the phase-space variable A(A) is given by /3(A)/3A. The equation that deter- 
mines the running of A follows from, (|3.5| ): 



A' = jX{X)(3 + 2& )i (18a + (2b 2 + 36 x ) log(l + A 2 )) ^ . (5.3) 

We give the plot of the potential in fig []]. The solution of ( |5.3| ) as a function of u, 
with the initial condition A(0) =0 is shown in fig ^. 

We see from the figures that A diverges in the IR at the singularity u$. 

The associated geometry can be derived as follows. The metric in the string 
frame is, 

ds 2 = e 2A+ ^dx 2 + e^dr 2 = \i(e 2A dx 2 + (A')~ 2 rfA 2 ). (5.4) 

The metric function A is obtained as a function of A from ( |4.18|) . Unfortunately, it is 
not possible to express this integral in terms of elementary functions. The Einstein 
frame metric is obtained from (|5.4| ) by multiplying with A~ 4//3 . The asymptotics of 
this geometry are discussed in the previous section. 
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Figure 2: The running of the coupling constant. It diverges at uq. 

One can also easily determine certain invariants of the geometry. The Ricci 
scalar is given by ( |3.9| ) and R = —3/2E. In the Einstein frame, one obtains, 

R = - 4X(A) 2 )(3 + 2b X)l (18a + (2b 2 + 3&0 log(l + A 2 ))^ . (5.5) 

The asymptotic behavior of the various geometric invariants in this geometry are 
also presented in the previous two sections. 

5.2 /^-function with an IR fixed point 

We choose an exact /3-function as 

1 + ^A - ^A 2 

with 

&i = &i - 26 2 /3 > 0. (5.7) 

In perturbation theory we have (3 = —b X 2 + 6iA 3 + C(A 4 ) with 6 01 > 0. 
According to (|3.8|) the potential for this specific /3-function is, 



= (A+)a+( l A )a _ (^+ - W - ^-) - ^(A - A) 2 j (A + -A) a +- 2 (A-A_) a - 

(5.8) 

Here we defined the following quantities. A = &0/&1 is the value of the coupling 
constant at the IR fixed point. A± is defined as: 



V 




Figure 3: Dilaton potential as a function of A when there is an IR fixed point. We have 
set Vq = bo = b± = 1 for simplicity. 



and the exponents are, 



4 At , 

0± = - T n-^:- (5 - 10) 



We note that along the RG flow, the range of A is such that, A_ < < A < A < A + . 

We give the plot of the potential in fig. |^. The RG flow is from A = in the UV 
toward A = Aq in the IR. As we discussed in section 4 the UV geometry is AdS with 



the radius given by eq. (|4.12| ). As fig. |3| suggests, the IR end of the potential is also 
AdS. 

Below we present the details of the IR geometry. The IR AdS radius is obtained 
by taking the A — > Xq limit in the potential, 



-IR 




V f = V ( ) . (5.11) 



We may deduce from (|5.11|) that Vf/Vo > 1. Therefore, the c-theorem of [36| is 
obeyed for this holographic RG flow. 

The equation that determines the running of A follows from, Q3.5| ) as, 

A' = ^|^A 2 (Ao - A)(A+ - \)^-\\ - \^-\ (5.12) 



The solution of ( 5.12 ) as a function of u, with the initial condition A(0) = is shown 
in fig |j. We observe that, A flows into a non-trivial fixed point at low energies. The 
value of the coupling constant at the fixed point is Aq = &o/&i- As A increases from 
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U 



Figure 4: Running of the coupling constant in the Banks-Zaks case. We have set Vq 
b = h = 1. 



to Ao during the RG flow, one has the inequality A_ < < A < Ao < A + in this 
range. 

The metric in the string frame ( |5.4j ) is obtained by using (|4.18|) and ( |3.5| ) as, 



ds 2 



2A - 



16 



dX 2 



^X b o (A -A) b o dx 2 + — (A -A)- i (A + -A)- a +(A-A_)- a -— . (5.13) 



V 



At 



One can also easily determine certain invariants of the geometry. The Ricci 
scalar is given by ( |3.9|) and R = —3/2E. In the Einstein frame we obtain, 



R = --V (5(A - A_) 2 (A + - A) 2 - 2A 2 (A - A) 2 ) (A - A_) a - (A+ - A) a + (5.14) 

Note that R goes to a constant in the IR, i.e. as A — > Ao- This is in accordance with 
the presence of the fixed point in the IR. Similarly the invariant (<9$) 2 follows from 
( |5.13|) as, 



Vc 



(9$) 2 = g x \d x ^Y = -^X\X - A_)-(A + - Ar+(A - A)' 

lb 



(5.15) 



which vanishes in the IR as required from the conformality. 

The IR asymptotics of the geometry can be obtained by solving the equations 
( |3.6|) and ( |3.5|) near A (that corresponds to u = oo). We present the geometry in 



the conformal coordinate, by performing the change of variable ( 2.15|) . One finds 

I 2 1 



(5.16) 
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near r — ► oo. Here A is an integration constant of equation (|3.5|) . The gauge theory 
interpretation of A is the energy scale below which the conformality in the IR sets 
in. 

We may determine the exact form of the quark potential and check that it indeed 
has the conformal form in the IR. We find the following relation between the quark- 
anti-quark potential and the distance in the IR asymptotics, 

E= C ^M\ (5.17) 

with the constant C is given by, 

1 5 

2h 2 A 3 

C[A °' bl] = TWn T^Tl x F (5 - 18) 

VdVo (A + - A ) 2 (A - A_) 2 

This exhibits conformal behavior. 

Of course, as the IR theory is conformal, we do not expect this model to be a 
good for description of pure YM type theories. We have presented it however, for 
illustration purposes. 

6. Fluctuations near the boundary 

The general geometry discussed in the previous section is asymptotically AdS near 
the boundary, u —>■ oo. It is important to study the fluctuations of the fields on 
this background. According to the standard rules of AdS/CFT, near the boundary, 
these fluctuations should correspond either to sources or expectation values of the 
corresponding operators of the gauge theory. In what follows, we compute the UV 
asymptotics of the fluctuations keeping the subleading terms as well. According to 
the principles of AdS/CFT these subleading terms should give information on the 
anomalous dimensions of the operators dual to the fluctuations. 

For the case of a 5D geometry coupled to one scalar field the analysis of the 



fluctuations have been carried out in |37|, |38| (see also [^] and f40|). We are interested 
in the fluctuations of the metric, the dilaton and the axion. The metric and the 
dilaton can be decomposed into 4D traceless-transverse tensor h^, two transverse 
gauge fields and and 5 scalar fields. 

The tensor modes obeys the usual 5D wave equation: 

□ 5 V = 0, (6.1) 

where CI5 is the scalar Laplacian on the 5D space-time 19 : 

0A ( d 2 dA d \ . 
V dr A dr dr J 



19 



In this section we use the conformal coordinate system, eq. (2.13) 
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From (|6.1|) , it is clear that the non-normalizable solution for h^ v sources a dimension 
4 operator. It is the energy-momentum tensor of the 4D boundary theory. Also it 
can be seen that the vector fluctuations can completely be gauged away by using 
5D diffeomorphisms. 20 Thus, in what follows we focus on the interesting case of the 
scalar fluctuations. 

6.1 Dilaton fluctuations 



As discussed in detail e.g. in [40|, there is a single diffeomorphism invariant combi- 



nation among the 5 scalar fluctuations: 



C = 1> - g^X- (6-3) 



Here X is the variable defined in ( |3.4| ), ip is related to the trace of the metric fluc- 
tuation and x is the fluctuation of the dilaton. We note from (|4.2| ) that, for any 
asymptotically-free boundary theory, X approaches zero in the UV. As a result, 
close to the boundary, ip and x decouple from each other and ( reduces to the dila- 
ton fluctuation. In the conformal coordinate system ( |2.13|) , it obeys the following 
equation: 

C+(3i + 2|)C-A; 2 C = 0, (6.4) 

where is the 4D momentum. 

It is interesting to comment on the correspondence of the gauge invariant scalar 
fluctuation (, with the Yang-Mills operators. Indeed, it corresponds to the "dressed" 
Tr[F 2 ] operator 

C <- (3(X)Tr[F 2 } (6.5) 

which is known to be RG invariant (see for example [0). The invariance of ( under 
infinitesimal r reparametrizations, guarantees that it corresponds to a RG-invariant 
operator. 

Using the asymptotic forms of the dilaton and the metric in ( |4.15| ) and ( [4.16! ) 



close to the boundary, (|6.4j ) becomes (in units where £ = 1), 



< - 1 (3 + i-L, + -g- - 2 t log '-;° E ; Ay ) C - ^ = 0, (6.6) 
r V logrA log 2 rA log 2 rA ' 



where b is defined in (|4.13|) and the coefficient C is given by, 

C = | + 26(2 - log 6o). (6.7) 
The corresponding fluctuation equation in the pure AdSs geometry is given by, 

C - -C - k 2 ( = 0, (6.8) 
r 



20 



This is strictly true only for the 4D massive modes, E0|, which is the case we are interested in. 
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There are two independent solutions to (|6.6| ), one is normalizable and the other is 
non-normalizable. The asymptotic expansion of the non-normalizable solution reads, 



(kr) 
4~ 
{kr) 4 
64 



1 - 



1 2-C 
log r A 2 log 2 r A 
3 



& log(-logrA) 
log 2 r A 



7-3C 361og(-logrA) 



2 log r A 4 log r A 



2 log r A 



+ C(r 6 X6.9) 



According to the standard rules of the AdS / CFT correspondence this solution is dual 
to a source term for a gauge-invariant operator of engineering dimension 4 and parity 
+ 1 in the boundary gauge theory. As argued above such an operator is /3(A) TrF 2 . 

The other solution of ( |6.6| ) is normalizable near the boundary and it has the 
following asymptotic expansion: 



(kr) 



log 2 rA + 2b log rA log(— log rA) 



26 - C - - ) l^/A - 



+ 0(r 6 ) 
(6.10) 



This solution corresponds to a VEV for TrF 2 , namely a gluon condensate. 



6.2 Axion fluctuations 

The other scalar in the geometry is the axion. Its background value is non-zero if 
Guv 7^ and is obtained by solving 



1 



— d» [Z($)y/ggVd v ] a 
V 9 







(6.11) 



in the Einstein frame. Z is an "effective" potential that resums higher-order correc- 
tions due to the five form, as explained in appendix |B.1|. The background solution 



a(r) describes the RG running of the QCD 9 angle. As the axion is suppressed by 
extra powers of 1/N C , f28|, this equation can be solved in the background already 



determined by the metric-dilaton equations. We analyze it in detail in |20| where the 
vacuum condensate of Tr[F A F] is calculated. Here we note that since the action 
for the axion is quadratic, its background value does not affect its fluctuations that 
are dual to the Tr[F A F] operator in QCD. 

The equation for the axion fluctuations follow from (|2.4|) and we obtain , 



a+ [3A + 2- 



k 2 a = 0. 



(6.12) 



The non-normalizable axion fluctuation corresponds to a source term for the 
operator TyF A F and the normalizable solution is dual to a VEV for the same 
operator, (TrF A F), namely an axial glueball condensate. This will be further 
determined in |2T| . 
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APPENDIX 



A. Asymptotic expansions in the (UV) weak-coupling regime 

In this appendix we will carefully analyze the weak coupling expansion in the UV 
regime of the space, near the boundary. We start from the equations of motion in 
the conformal frame as described in section 



ds 2 = e 2A(r) [dr 2 + rjijdx^x^ , E„ = QA 2 , E {j = 3 



A + A 



Vij (A.l) 



where we use the mostly plus convention for the metric. The two scalar invariants 
of the metric up to quadratic order are 



E M 



-2 A 



2A + 3A< 



, E^E^ = 36 e~ 4A 



A 2 + 2AA 2 + 2A 4 



In this ansatz the equations ( |2.11| ) become 
4 • 



4 „ 



12A 2 - -0 2 - e 2A V = , 6A + QA 2 + -0 2 - e 2A V = 
3 3 



(A.2) 



(A.3) 



+ 3A<p + 3 -e 2A V> = , V > = ^± 



(A.4) 



To continue further we define the 't Hooft coupling and its inverse as usual 



A = N r e 4 



A 



o:.s 



a, at 



(A.5) 



( [OP and (|Al|) become 



a, at 



• a 



3A— + ~-e 2A Y = 



3 a; 



e 2A ^ = 



(A.6) 



We now change variables to 



r = e 



e A ^ = e'd(t) 



(A.7) 



where t is essentially the log of the energy in the perturbative region. The equations 
become 



- — - 4^ + -4- - 3—— + -d 2 V = 



a. 



12 1 



1 



a. d 8 



3 aj 



(A.8) 
(A.9) 



where we now use primes for t derivatives. 
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We assume a regular power expansion for the potential 

V = V + V 1 X + V 2 X 2 + V 3 X 3 + V A X 4 + --- , V = V 1 X + 2V2X 2 + 3V 3 X 3 + 4V 4 X A + ■■■ 
(A.10) 

We now substitute for the inverse coupling 



1 

x =as 



h foflogL (b\ b 2 \ 1 6 3 log 2 L 

L -V^ gL + l?— + [^ + VjL + 2bJ—- 



(A.ll) 



bib 2 log £ 
6^ L 2 + 



26 26g 



1 ^ / 1 



L = a + b t = -b log(rA) , L' = b 



(A.12) 



The expression above is taylored so that coupling satisfies the standard RG equation 



A' = -b X 2 + M 3 + b 2 X 4 + 6 3 A £ 



, h b 2 b 3 
b n , + 



a q a? a? 



(A.13) 



where A is the RG invariant scale of QCD. 

We may rewrite now equations ( |A.8| ), ( |A.9| ) as 



d! la s ^ 2 y/ 4 la" 1 a' s d! 
d 8 a' 3 3 a' 3a s d 



<PV 1 a? 



(A.14) 



so that 



i^dV - - - + 

8 a' 3 3 «'„ 3 a. 



W laf 
~12~ + 9 a? 



(A.15) 



Note that only one branch of the square root is relevant. We expand the potential 
in inverse powers of logs to obtain 21 



V x V 2 6i logL / 6 2 6? \ 1 



(A.16) 



a's b \b bl J L 
From ( |A.15 ) we obtain 



bo bl ) L 2 ' 



2b,V 2 bJVA 1 + logL (l_ 

bi + bi ) i? +u U 3 

(A.17) 



3 -d 2 



V + 2*V' (1 + ^-2* 



a' a. 



(A.18) 



21 We thank Shesansu Pal for pointing out errors in the following two equations in the first version 
of this paper. 
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V + <*fV 1 + 2f- - 2* 
1 a' a. 



a' / \ a 



( ZfV 



(A.19) 



Using 

a 



Srrf 



ad 



60 



+ 



3V 3 , 26iV 2 

"7 1 Ta v<2 



(2V 2 




V 6 


bl 


b 2 V 1 
b 2 


biVA 
bl J 



+ 



b 2 



1 



+ 



b\ ) L 2 



(A.20) 



we obtain 



boV \ (2b Q V + V 1 )(b 1 V 1 + 2b V 2 ) b 1 {2b Q V + V 1 ){b 1 V 1 + 2b V 2 )logL 



Vi 



boVfL 2 



+ 



(A.21) 



, 62(260^ + ^) , 61(61^ + w 



V?L 2 



± 



V 2 L 2 



V 2 L 2 



± 



(2b b 1 V 2 V 1 + ShVoV 2 - boV* + Ab 2 V 2 V 2 + 660^0 W 2 ) , ( \ 



^ 2 V' U V W v ~» V^V (2V l + b V )L 

Compatibility to leading order with the first of ( |A.14j ) implies 



+ OI- 



Vi = lb V 



(A.22) 



for the plus sign. The minus sign does not lead to a consistent solution. Therefore 
to leading order 

2 



3 J2 V _ b o y o 



1 + 



60^0 \ ( b V Q 



Using ( |Q2D , ( [70TD becomes 



1 + 2-% = ^ =► dV = 12 (A.23) 

Oo^O * 



9 3(6^-126^0-27^) 36 1 (-6 V + 126^ + 27V 2 ) 



(A.24) 



2 56 VoL 56 2 K)L 2 

(646q- - 75 66 2 6x - 19446 2 - 54 006 6 2 )\/ 2 + 23496^0^ - 208986^0^2 - 37179Tg + 182256oVftV3 

75062\/ 2 L 2 



+0 



1 

U 



We may calculate from ( [A.24 ) 



d! b 2 V - 126i V Q - 27V 2 2b 1 (b 2 V -12b 1 V -27V 2 )\ogL 



d 



IW L 2 



156 K)L 3 



(A.25) 
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^946 4 - 12516^6i - 3246? - 54006 6 2 9(36g - 316i)V 2 



567U 9 2 



33756 £ 3 125& V L 3 
Compatibility with the first of (|A. 14|) implies 



27V 3 



125b v 2 L 3 5V L 3 



V 2 _ 236§ - 366i 
V ~ ¥ 



Vs 3246 2 + 124&g + 189&0&! 
V Q = 37 



(A.26) 



V 4 _ 31766^ + 72366§6i - 2436? - 8106 &2 - 43746 3 
V ~ 



3 9 



Using this we finally obtain for the metric 

72 QL A(r,RU2 



12 " £ 2 



86 4(266^ + 9h) 



86i log L 8 b\ log 2 L 



3 4 L 2 



3 2 L 2 



3 2 6 n £ 3 



16&i(136g + 9&i)logL 8 



u :f b i 698bi o + 594& o&i - 2436 2 - 3246 6 2 )— + O 
Therefore to next to leading order the Poicare metric near the boundary is 



(A.27) 



1 

I 1 



ds 1 



1 + 



3 2 logrA 



+ 



4(26 + 9^-18fHog(6 log-±r 



3 4 log 2 r A 



+ o 



log logrA 
log 3 r A 



- (dr 2 +dx 2 ) 
'(A.28) 



We observe that the first non-trivial correction is independent of bo- This is related 
to the fact that as we have no unambiguous identification of the normalization of the 
gauge coupling, the /3-function coefficients we are using can be changed as 



(A.29) 



From (|A.27|) we observe that d is invariant up to the shift of the logarithms, log L = 
log(&olog^-) since it depends only on b\jb\ and b 2 /bl- It is a non-trivial statement 
that a change of scale in the coupling as in ( |A.29|) , in ( |A.27| ) can be absorbed into a 
change of the scale A as well as the non-universal /3-function coefficients b n> \. 

The potential in the UV regime can be therefore expanded in terms of the overall 
AdS scale as well as the /3-function coefficients to obtain 



V 



12 



1 + g(M) + 



23 - 36^ 



3 1 



HboXY 



324|§ + 124 + 189^ 



3 7 



^(M) 3 + ^(A 4 ) 

(a:3o) 



A.l Scalar curvature invariants in the UV regime 

We may now evaluate the basic scalar invariants of the background in the weak 
coupling (UV) regime: 



V = Vn 



8b 23b 2 861 log L 

1 H H H — h O 

9 L 81L 2 9L 2 



1 

Z 3 



(A.31) 
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m 2 = e~ 2A % 



■ 2 



r 2 d s 2 
d 2 a 2 



1 a' 2 

s _ y 

d 2 a 2 



b | 6 logL | 86§-96i 



12L 



12L 2 



108L 2 



R = _(<90) 2 - -V 



(A.32) 
(A.33) 



3 N ~ T ^ 3 1 

It is obvious that all invariants are regular in the UV, ( L — ► oo) . All higher curvature 
invariants are also regular as they are polynomial functions of 



j(8*)" + ^ , e -»A=->)' + l V 



-2A^2 _ 1 /qj.\2 ' 



(A.34) 



A. 2 The two-loop /3-function coefficients of pure gauge theory 



As a final point we quote the scheme-independent /3-function coefficients for QCD, 
b ,h. 

With Nf (non-chiral) flavors in the fundamental, the /3-function reads 



(4tt) 



-a, - In* > 



(4vr) z 



2*a? - *> 

3 c AT„ 



For the the 't Hooft coupling, and with ^ 



3 " c 
a; we obtain 



+ 



(A.35) 



X = g 2 N c , A 
from where we obtain 

b = 



ll-2x) x2 , (34-13x) x3 



(4vr) 2 

2 (11 - 2x) 

3 (4vr) 2 ' 



In this paper, x = and therefore 

°0 



61 

1.2 
°0 



3-34 
2-121 



(4tt)« A + 

3 (34 - 13x) 
~2 (11 - 2a;) 2 

0.42. 



(A.36) 



(A.37) 



B. General potentials at string tree level 

Extra terms in the dilaton potential can arise from higher a' corrections proportional 
to the field strength of the RR four-form. We parametrize these corrections in the 
string frame as 

Sf = ~^ J d 5 x^ e^K (e^y) , (B.l) 

where 

^— ' n 5! 

n=l 

and we assumed the simplest type of contraction for simplicity. The contribution for 
general type of contractions yield qualitatively the same result upon modifying the 
coefficients, a n . Therefore we will assume that ( |B.1| ) captures all such corrections. 
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Passing to the Einstein frame we obtain 

M 3 



Op 



2tl 



J d 5 x^/ge^K (eT^y 



(B-3) 



The equations of motion are 



7/«l 



e-^K' ( e -*M W-m. 



(B.4) 



with solution 



= Ne**E 



(B.5) 



where N c will now be non-linearly related to the number of color branes. By squaring 
we obtain the consistency conditions 



y 



K'(e-%+y 



-N 2 c ef* , (K') 



l\2 ^ 



N? 20 



4! 



<t> n 



(B.6) 



so that 



(B.7) 



We now compute the contribution of the five-form to the stress tensor 



TF = ±SSl 



M 3 



4! 2£? e 



M 3 ei< 



y 2 



We may therefore substitute the action of the five form Sf with 

M 3 f 



V E 



eh 



In the formulae above y = X(4>) is a solution of ( |B.6| ). 
We redefine 

C = e"-^ 

Then the relevant equations become 



9iiv 
(B.8) 



(B.9) 



(B.10) 



fB.ll} 



C[K'(C)] 2 = -N 2 c e 2 * , V E = e 



CK'(C) - \k{0 



(B.12) 



We return to the a-model frame in order to estimate the large- N c dependence 



S a = M 3 J tfxyfge 



2<p 



i? + 4(<90) 2 + -K(0) 



(K>(() - -K(C) 



(B-13) 
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If we define the 't Hooft coupling as 



A = N c (B.14) 

then from ( |B.12j ), ( is a function of A. The a- model frame action ( |B.13[ ) becomes 

(d\) 2 1 



5. 



(B.15) 



We observe the 0(N 2 ) dependence expected from the sphere action as well as the 
fact that there is a non-trivial potential for the 't Hooft coupling. 
For small ( we can calculate the potential perturbatively: 

K(() = -2Sc + ( + ±a 2 ( 2 + 0(( 3 ) , (K'(0 = ( + a 2 e + O(( 3 ) (B.16) 
to obtain 

C = -A 2 - 2a 2 A 4 + 0(A 6 ) , V a = -\ (-25c + A 2 + ^A 4 + 0(A 6 )) (B.17) 



For a DBI-like example : K(y) = 2y/l + y we obtain 



A 



2 



C = - TT ^ , V(X) = -VT+>? (B.18) 

It is obvious from the above analysis that V(A) has an infinite series of terms, 
that, after solving the equations of motion of the four-form, are independent from 
a' except an overall (universal) dependence in front. Therefore the higher-order in 
a' corrections to the four-form are equivalent to a (leading) order in a' potential. 
Surprisingly, this potential has a regular expansion at weak YM coupling. 

This analysis obviously generalizes to the kinetic terms of the dilaton and gravi- 
ton. The general action is of the form 

S = J d b x^e- 2 ^ [V R+ ^Vi{d(j)) 2 + V 2 ] (B.19) 

The functions V, V\^ 2 depend on e 2 ^y with y given in ( |B.2|) , and summarize the higher 
a'-corrections of the four-form. Going through the same procedure as above, it can 
be shown that this is equivalent with integrating out the four-form, and multiplying 
the kinetic terms with appropriate series in the exponential of the dilaton as in the 
standard potential. Moreover, higher derivative terms proportional to powers of the 
curvature as well as the dilaton derivatives will be generated. 

B.l Corrections to the axion terms 

The axion a, dual to the instanton density has a special position among the super- 
gravity fields. Its shift symmetry is protected, and its special normalization in the 
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YM theory implies that its contributions are suppressed by a power of 1/iV 2 , due 
to the fact that the dual variable in YM is an angle, [p8| |. This is reflected in the 
fact that in string theory, a is a RR field and therefore has suppressed dilaton de- 
pendence. The leading term in the effective action is J ^J~g (da) 2 in the string frame, 
and becomes 

J Jge 2 \da) 2 = ± J ^\ 2 (da) 2 (B.20) 

in the Einstein frame where in the second equality we indicated the suppression of 
the ^-induced vacuum energy. 

We now consider the higher terms in the a 1 expansion that involve the 5-form 
field strength. In analogy with the previous section we may write them in the string 
frame as 22 



M 3 
2£2 



d 5 x^ e" 2 * [K x (e^y) + K 2 (e 2 *y) e 2 *(da) 2 ] 



(B.21) 



where 



Kl (y) = -26cJ2^V n , K 2 (y) = l + J2 



71=1 



n=l 



J n n 

—y , y 

n 



5! 



a i 



Passing to the Einstein frame we obtain, 

M 3 



qE _ 

211 



dfxy/ge** Kx le"^ + K 2 (e-^+y) e^(da) 



The equations of motion for the 5-form are 



e-f+K[ (e-^y) +e -^K 2 e^M (da 



with solution 

K (e-%+y) + e&K' 2 (e~^y) (daf 
By squaring we obtain the consistency conditions 



F =0 



F = N e^^E 



y 



K'x (e--*y ) + e^K' 2 ( e"M (da 



-N 2 e~ 



and 



so that 



K[ + e^K 2 (dc 



1 2 F, 



4! 

F 2 



N 2 C 20, 



flU 



(B.22) 
(B.23) 

(B.24) 
(B.25) 

(B.26) 
(B.27) 
(B.28) 



22 Higher terms like (da) 4 are suppressed by extra powers of N c and we do not need to consider 
them here. 
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We now compute the contribution of the five-form to the stress tensor 
2 „„ 



M 3 ' 



2 1 6S§ 

- e s < 



2if 



If" 



14 ± 

e-^y 



K[ + e*+K' % (dc 



Kx + e&Kzidaf 
K x + e^K 2 (da) 2 



Fl, e 2 * 



|; + -^-K 2 d /J ,ad v a- 



(B.29) 



=20 



The axion equation also reads 



We redefine 

and ( |B.26| ) becomes 



V [e 2 ^ 2 V M a] = 



crac)+^ a (Ch] 2 = -^ 

The dual action which gives the same equations of motion is 



fiW + -^K 2 d ll ad v a 
(B.30) 

(B.31) 

(B.32) 

(B.33) 



~S E = ^- I <?xy/ge 



am) + v k' 2 (0) - -(^(o + ^ 2 (o) 



(B.34) 



In the action above, C{X,f]) is a solution of ( |B.33[ ). In the string frame the dual 
action becomes 



((^(O+raj-i^W + M)) 



(B.35) 



with 77 = e 2 ^(da) 2 . We must separate the kinetic term of the axion from the higher 
derivative terms that appear because ( depends non-trivially on rj. This turns out 
to be in the Einstein frame 

M 3 



on, 



axion— linear 



u 2 



d 5 x^ e 2 <f>K 2 (C)(dc 



(B.36) 



where now (* is a solution C^i(C) = —N^e 2 ^. 

We finish this section by giving the most general dual action involving the five 
form 

S = J d 5 x^/g e- 2 * Z(e 2 *y, Zi ) (B.37) 

where Z is an arbitrary function and z% are scalar invariants of other fields, R, R 2 , R tlv R tJ,v , 
(<90) 2 , e 2 *(da) 2 ) etc. 

The dual action is given by the Legendre transform 



(d ( Z((,z t )-^Z((,z t 



where ( satisfies 



C(d ( z((,z l )) 2 = -\ 2 



(B.38) 
(B.39) 
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B.2 Corrections to the gauge coupling constant identification 

Consider now a probe D3 brane and the coupling of the kinetic gauge field terms 

Sd 3 = ^ / d^ge-*Z{e^y)Tr[F% Z(y) = 1 + £ y = 3^ (B.40) 

s n=l 

where T3 is dimensionless and Z(y) summarizes higher-order couplings of the five 
form on the D-brane (arising from disk diagrams). Going to the Einstein frame and 



substituting from |B.ll| , pTT2| we obtain for the gauge coupling constant 

9ym = (B.41) 

Therefore, these corrections, although without derivatives, are due to the higher a' 
corrections on the branes. A similar argument indicates that the tension of possible 
flavor branes obtains similar types of corrections. 

B.3 Other higher derivative corrections at the tree level 

After dualizing the five form, there remain true higher-derivative corrections asso- 
ciated to curvatures and derivatives of the dilaton. We will indicate here how such 
corrections can be instrumental in effectively generating a constant term in the dila- 
ton potential, a fact that we have assumed in this paper in order to simplify our 
problem. 

To give the idea, we will focus on the tree level, string frame effective action with- 
out dilatonic kinetic terms and considered only as a function of the scalar curvature 
(for simplicity) 

S = J d'x^e^ f(f s R, e 2 ">y) (B.42) 

where we are using the same definition of y as in (|B.2|) and / is a function with a 
regular series expansion around R = y = . The equations of motion are 

i 2 sf R R,u-lg,uf + e^[g, u D-V,V u ](e^f) = , fn = ^ (B.43) 

v"[/»W = o > f y = jjj (R44) 

These equations admit an AdS§ solution, 

4 

R v» = -12—9^ (B.45) 

l AdS 

with =constant and 

f yFvv .. v , = ^E Vv .. VB - y ft [x, e 2 *y] = -JVj (B.46) 
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where we defined 

x = -2o| (B.47) 

The metric equation then becomes 

2xf R [x, e 2 *y] + 5/ [x, e 2 *y] = (B.48) 

and typically the two algebraic equations are expected to have a solution. A simple 
example involves 

which generates an AdS$ space with 

e 9/i 



P s 10(7 ± v/49 + 180/xA 2 ) 
In particular, a solution exists for A = and it is 

I 2 _ 9^ 

e ~ i4o 



A = e*N e (B.50) 



(B.51) 



It is not difficult to verify that allowing now the dilaton to run logarithmically 
in the UV, so that A^Oasr^Oisa small (subleading) perturbation in the a' 
expansion, generating a solution that logarithmically asymptotes to the AdS solution 
above. In particular, the kinetic terms of the dilaton that have been neglected will 
be suppressed by large inverse logs, and the same applies to the power series in the 
t' Hooft coupling N c e^. We will not pursue this avenue further here, leaving it for a 
future investigation. 

C. Perturbative analysis near an extremal (AdS) point of the 
dilaton potential 

In this appendix we will analyze here the stability properties of the perturbative 
dilaton potential. 

We will use the domain- wall coordinate, ( |2.14| ). The field equations are given by 
( |2.18| ), ( [2.18] ) and (|2.19| ). We parametrize the potential around an AdS extremum as 



where $ << 1. Since V f = at the critical point, there is a AdS§ solution with 

A=j , $ = (C.2) 

Perturbing around the fixed point solution A = A* + 8 A, $ = 5$ we obtain to linear 
order 

^■5 A' = <5$ /2 - J-$ 2 = 0(<5$ 2 ) , 5&'-^5&-^5$ = (C.3) 
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Figure 5: Left:The dilaton potential plotted as a function of the t 'Hooft coupling for 
x = 5. Right: The value of the t 'Hooft coupling at the extremum as a function of x. 



We observe that to linear order 5 A is a constant, which amounts to a renormalization 
of the AdS length scale I. We can therefore ignore it. The general solution of the 
second equation is 

5$ = C + e t + C_ e t (C.4) 
Changing variable to r = Le u ^ we obtain to linear order in the perturbation 



ds 1 



" 2 t\ ( 2 +Vl+?) /r\( 2 - 2 Vi+?) 



; (dr A + dx 2 ) + ■ 



5$ = C 4 



+ C_ 



(C.5) 



We may now analyze the potential (we assume we are in 5 dimensions and 
therefore 5c = 5) 

.1 

1 



5 - -A 2 - x X 
2 



where 



x = 



N c 

The potential is plotted in figure || 

It has a single extremum at A = Aq with: 



An = N r e 4 



-7x + V49X 2 + 400 
10 ' 



From this we obtain the value of the potential at the extremum that gives, 



P2 \3 
-s A 



400 



100 + 7x 2 - xV49x 2 + 400 
400 



as well as the second derivative parametrized as in (|C.1J) 

5 



4 



400 + 49x 2 - 7xVA9x 2 + 400 



100 + 7a; 2 - x^/A9x 2 + 400 



(C.6) 
(C.7) 



(C. 



(C.9) 



(CIO) 
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L 




Figure 6: Left:£ plotted as a function of x. Right: The value of the AdS radius L in units 
of t s as a function of x. 

As x varies between < x < oo, the parameter £ that controls the anomalous 
dimension of the YM coupling constant varies as 5 > £ > | as can be seen in figure 

B 

This implies that the associated dimension is A = 2 + 2y/l + £ and satisfies 
2 + 3a/2 < A < 2 + 2 a/6 or equivalently 6.24 < A < 6.90. It therefore corresponds 
to an irrelevant operator. This is not what we expect for QCD at weak coupling. 

The AdS radius i (in unite of l s ) becomes small at x < 0, whereas it increases 
as x becomes large as can be seen in figure [| Therefore it becomes arbitrarily larger 
than the string scale for large x and therefore the background will be weakly curved. 
Therefore, although the t 'Hooft coupling at the AdS extremum is small as x becomes 
large, the geometry is weakly curved. This is against the AdSsxS 5 intuition. 

It should also be noted that the coefficients N c and Nf appearing in the potential 
are multiplicatively related to the number of colors and flavors respectively. They 
equal them if the relations stemming from M = 4 branes still hold, but this is not in 
general guaranteed. 

D. Analysis of the dilaton potential with single exponential 

In the string inspired dilaton potential, the weak coupling asymptotics of the flow is 
governed by a single exponential in the potential. Also in the confining examples that 
we study in this paper, the leading term of the potential in the IR is an exponential. 
Therefore it is desirable to investigate the solutions of the system given by a potential 
of the following form (after a convenient shift in the dilaton) 

V(<f>) = ^ee a *, 6 = ±1. (D.l) 

We want to find all the solutions of the system of equations, ( |2.18| ) and ( p.lSQ 
with V given by ( p.l| ). As a starting point, we classify the solutions by the behavior 
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of the phase space variable X which, in this case, obeys the following simple equation 
obtained from ( |3.7|) : 

« + (D. 2) 

d(f) 3 a 6X 

where we defined 

a = 8/3a. (D.3) 

The fixed points of ( [L).2| ) are given by X = ±1 and X = — 1/a. 
Equation (|D.2|) can be integrated to yield 



(X- + l)w=v 



3a 

2T 



(X + ±)«a--' 
For the special case of a = ±1 we obtain instead 



;d.4) 



log fTT-xTT = f w - w - (a5 > 

Now, the solutions of the system are given by the fixed point solutions of (|D.2| ) 
and the solutions that flow between these fixed points. We first list the fixed point 
solutions: 

D.l The X = 1 fixed point 

This amounts to solving 

<f/ = 3A' , 0" + l0 /2 = O , e a(p ^0 (D.6) 



with solution 



e* = C(u -uf* (D.7) 



We take the range of u as u G (— oo,u ). If a > this solution is valid in the 

immediate neighborhood of u = uq. If a < it is valid in the neighborhood u — > oo. 
We also find 

e A = C(u - u)* (D.8) 

and for the Poincare coordinate 

i 

r-r = ^(« -w) 3 , & W=(^j (r-ro)s , = -C(7(r - r ) (D.9) 

This solution is valid near r = r Q . 

D.2 The X = -1 fixed point 

This amounts to solving 

<f/ = -3A' , <p" -^0 ,2 = O , e a *-^0 (D.10) 
It is related to the X = 1 solution by </> — >• —(f). 
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D.3 The X = -\ fixed point 

It exists if 

(a) e = 1 and \a\ > 1 

(b) e = — 1 and \a\ < 1 

(c) e = ±1 and \a\ = 1. In this case it merges with the X = ±1 fixed point 
solutions. 

The solution is 

_ 3a 

e*=(— — ) * , e A = C(u -u)^ , C=^v/ e (a 2 -l) (D.ll) 
In the Poincare coordinate, 



2 



r ~ r ° = (4-o 2 )C ~ ' = ' " ' ' '' ~ ' ° ' "' 



3a 

e* = C-^l {A -f )C Y~ a \r-r )^ (D.13) 
D.4 Flow solutions 

More general solutions are given by the flows between the fixed points of X. The 
nature of these flows are determined by the stability properties of the fixed points that 
we depict in fig.([7|) for the case a < 0. A cross denotes an unstable fixed point and 
a point denotes a stable one. The arrows between these points depict the direction 
of the flow in an appropriate radial variable that we define below. The boundaries 
of the phase space are given by X = ±oo and X = 0. In particular the X < and 
the X > solutions are disconnected. The case a > can be obtained from fig.(0) 
by utilizing the symmetry of (p.2|) under X — > —X, <p — > —0 and a — > —a. 
The system can be solved analytically in the following radial variable, t: 

4- = e^4- (D.14) 
du dt 

We describe the solution for the case e = +1 which is the interesting case for 
us. In this case the range of the phase space variable is given by — 1 < X < 1. This 
range is divided into two flow parts by the presence of the fixed point at X — — 1/a. 
Let us take a > 1 (a < 1 is given by a simple replacement of the hyperbolic functions 
below with the triangular ones). 

D.4.1 The solution: -1 < X < -1/a 
Let us define, 

t = ^- (t-t ) D.15 

3a 



(D.12) 
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Figure 7: The stability properties of the fixed points. Case A is for < a < 1 and case 
B corresponds to a > 1. The cases where a < can be obtained from the above figures by 
using the reflection symmetry in X. 

The solution is given by, 

3ft f 1 X "\ 

= 0o + -H — -r log cosh(t) logsinh(t) , (D.16) 

4 \a + 1 a — 1 / 

and 

A = A + I ^—j-y log cosh (t) + _L-logsinh(i)^ . (D.17) 

The phase space variable as a function of t is given by, 

^tanh 2 (t) - 1 

X{t) = ^ . (D.18) 

^tanh 2 (t) + l 

In these equations, we choose t > that corresponds to the range t < t < 00. 
This solution corresponds to the flow from the fixed point X = — 1 at t = t to the 
fixed point X = — 1/a at t = 00. 

The asymptotics of this solution are as follows. As t — > t , (X — > —1), 

A^O-to) 3 ^ (D.19) 
ds 2 -> (t-io)*^da; 2 + (t-to)^<ft 2 (D.20) 

Thus the space shrinks to a point as one approaches to the X = — 1 fixed point. As 
t -> 00, (X -> -1/a), 

A -> v^T (D.21) 
ds 2 -> — (dx 2 + r"^dr 2 ) (D.22) 
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where in the last equation we changed the variable as r = exp(— atj\Ja 2 — 1), hence 
r — > in the limit. We observe from the asymptotics of the metric, that the space 
becomes AdS only when a is taken to oo. This nicely parallels the fact that as 
a — > oo the fixed point approaches to A — > and, as our primary example in this 
paper suggests, one has an asymptotically AdS space in the A — > limit. 

D.4.2 The solution: -1/a < X < +1 
The solution is given by, 

<f> = 0o + x ( ~~T !ogcosh(t) + — ^- logsinh(t) ) , (D.23) 

4 \ a — 1 a + 1 / 

and 

A = A + - (—}— logcosh(t) + — !— logsinh(t) ] . (D.24) 
4 \a — 1 a + 1 / 

The phase space variable as a function of t is given by, 

-f±±tanh 2 (i) + 1 
2±itanh 2 (I) + l 



1(1) = -£=± ^ . (D.25) 



This solution is valid in the whole range — 1/a < X < +1 (for a > 1) and describes 
the two flow solutions: one from X = to X = —1/a and from X = —1/a and one 
from X = to X = 1. 

The interesting asymptotics are i — > oo, (A — > —1/a) and £ — > t\ where, 



tx = t + — 3a tanh^W- -. (D.26) 

2 v / o^T Va + 1 v ; 

In the first A — > —1/a, one obtains the same asymptotics as in ( |D.21| ) 

and ( p.22| ). In the second CcLSG, cLS A — > 0, one finds, 

A -> const. (D.27) 
ds 2 -> dx 2 + dt 2 . (D.28) 



E. function with an exponential tail 

In this appendix we provide another example of an exact /3-function and the associ- 
ated holographic geometry. We choose the /^-function here as: 

p = -b X 2 e ~ CoX (E.l) 

In perturbation theory this describes an asymptotically free theory. In the IR it 
reaches a fixed point but in way which is different from the previous example. In 
particular the fixed point is at infinite coupling. 
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Figure 8: Scalar potential as a function of A for the exponential type running. We set 
Vo = bo = cq = 1 for illustrative purposes. 



The potential that follows from this /^-function is obtained from (|3.8| ) as, 

V = V (l - fAV^e-S^- 1 ), (E.2) 
9 

The plot of the potential is shown in fig. [8]. The RG flow is from the Gaussian 
fixed point in the UV to an IR fixed point at A = oo in the IR. The corresponding 
UV geometry is an asymptotically AdS space with radius given by (|4.12) as before. 



The corresponding far IR geometry is also AdS as expected from the fact that V in 
(|E.2|) approaches to a constant as A — > oo: 



V -> V f = V e^. (E.3) 
The radius of the AdS in the IR is given by, 



The running of the 't Hooft coupling is determined by the following differential 
equation that follows from 



A' = ^ '-boX'e-^e-—^ L K (E.5) 

6 

We display the numerical solution to this equation in figure |^. From this figure it 
is clear that the coupling constant diverges in the IR. One can calculate the scalar 
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u 



Figure 9: Running of the coupling constant in the exponential case. We have set Vq 
b = c = 1. 



invariants of the geometry. The Ricci scalar is given by (|3.9| ) and R = —3/2E. In 
the Einstein frame, one obtains the IR limit value of R as, 

R = -^V f , (E.6) 

where Vf is given by ( |K.^ ). Note that the function X in ( |3.4j ) limits to zero both 
in the UV and the IR. This is in accord with the fact that the scalar potential goes 
over to constant values both in the UV and the IR. The 5D geometry is of the form 
of a Minkowski domain wall that interpolates between two AdS geometries. 

The details of the approach towards the IR AdS geometry can be studied as 



in section |5.2| . One determines the following subleading behavior near the IR fixed 
point: 

ds 2 = _ 4 1 1 } _ 2 , d 2 d 2 X ( } 

r 2 v 9rAloglog(rA) ; v ; ' y ' 

as r — > oo. 
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